00 

(V 


USMMRDL  TECHNICAL  REPORT  72-2 


rH 


G5 

<d 


INVESTIGATION  OF  THE  VORTEX  NOISE 
PRODUCED  BY  A  HELICOPTER  ROTOR 


By 

N.  Kevii  JihRSon 
Walter  M.  Katz 

Felriary  1972 


u.  s.  ARMY 


i 


EIISTIS  DIRECTORATE 

AIR  MOBILITY  RESEARCH  AND  DEVELOPMENT  LABORATORY 
FORT  EIISTIS,  VIR6INIA 


CONTRACT  DAAJ02-70-C-0023 
ROCHESTER  APPLIED  SCIENCE  ASSOCIATES,  INC. 
ROCHESTER,  NEW  YORK 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

Springfield,  Va  27151 


Approved  for  public  release; 
distribution  unlimited. 


DISCLAIMER  NOTIC 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


DISCLAIMERS 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the 
Army  position  unless  so  designated  by  other  authorized  documents. 

When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data  is  not  to  be  regarded 
by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person 
or  corporation,  or  conveying  any  rights  or  permission,  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 

Trade  names  cited  in  this  report  do  not  constitute  an  official  endorsement  or  approval 
of  the  use  of  such  commercial  hardware  or  software. 


DISPOSITION  INSTRUCTIONS 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to  the  originator. 


MCIBIM  far 

CftTI  WHITE  SrCTlOH 

IOC  IUFT  SECTION  □] 

IMKXOlMCn  □ 

jibtihmtio«  . .  i 


IT  ^ . 

IISTlltyTIM'mnASILlft  OMR 

list.  !  SMIL  Ul  «  VEUM. 


Helicopter  Noise 
Acoustics 
Fourier  Transform 
Vortex  Noise 


Filter  Bandwidth 


Unclassified 

Security  Cli«*iNc«tloa 


Unclassified 


Security  CIihUIciIIw 


[  DOCUMENT  CONTROL  DATA  -RID  1 

|  fSicuriff  ot  tlllm,  body  ot  obatrocl  mnd  Indaalng  annotation  mutt  bo  ontormd  whon  tbo  ommII  ropott  It  e  lotalilod)  | 

1.  ORIGINATING  ACTIVITY  (CotpOftO  OUtbof) 

Rochester  Applied  Science  Associates,  Inc. 

140  Allens  Creek  Road 

Rochester,  New  York 

U.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

J6.  CROUP 

1.  REPORT  TITLE 

INVESTIGATION  OF  THE  VORTEX  NOISE  PRODUCED  BY  A  HELICOPTER  ROTOR 

«.  OESCRIPTi  VG  notes  (Typo  ot  ropott  ond  Inctuahro  dm  ft) 

Final  Technical  Report 

•  ■  autmorisi  (Flrat  naata,  mlddtm  Initial,  taal  naata) 

H.  Kevin  Johnson 

Walter  M.  Katz 

•  REPORT  DATE 

February  1972 

70.  TOTAL  NO  OP  PAGES  76.  NO-  OF  REFS 

132  18 

M.  CONTRACT  OR  OR  AN  T  NO. 

DAAJ02-70-C- 0023 

6.  PROJEC  T  NO. 

Task  1F162204AA4104 

c. 

da 

OO.  ORIGINATOR*!  REPORT  NUMVE  R|SI 

USAAMRDL  Technical  Report  72-2 

S6.  OTHER  REPORT  NOIS)  (Any  otfiR  nuaib«f«  tfi«f  mmy  bo  oaalgnod 
thlo  roport) 

RASA  REPORT  71-10 

10.  DISTRIBUTION  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

II.  SUPPLEMENTARY  NOTES 

13.  IPONIORINS  MILITARY  ACTIVITY 

Eustis  Directorate,  U.  S.  Army  Air 
Mobility  Research  &  Development 
Laboratory,  Ft.  Eustis,  Virginia 

Karman-street-type  vortex  shedding  from  a  lifting  surface  was  analyzed 
as  a  source  of  noise  from  a  helicopter  rotor  in  hover  and  forward  flight. 
Experimental  pressure-time  histories  were  analyzed,  and  high  resolution 
spectra  were  developed  over  a  frequency  range  of  0  to  5000  Hz  using  a 
0.7-Hz  filter.  On  the  basis  of  these  spectra,  the  main  and  *:ail  rotor 
rotational  noise  and  discrete  noise  sources  were  identified  and  then 
removed  from  the  measured  pressure-time  histories,  leaving  pressure¬ 
time  histories  representing  only  the  broadband  noise  radiated  from  the 
rotor  system.  The  theoretical  analysis  that  was  developed  related  the 
acoustic  radiation  of  the  vortex  shedding  forces  on  the  blade  to  the 
experimental  broadband  pressure-time  histories.  This  analysis  was  then 
used  to  solve  for  the  oscillatory  forces  on  the  blade  that  would  dupli¬ 
cate  the  broadband  noise  characteristics  of  measured  pressure-time 
histories  of  the  helicopter.  The  results  of  the  investigation  indicated 
that  "vortex  noise"  is  the  major  source  of  acoustic  radiation  from  a 
helicopter  rotor  in  hover  or  low-speed  flight  and  that  it  is  concen¬ 
trated  in  the  frequency  range  of  200  to  500  Hz.  Because  of  the  excel¬ 
lent  correlation  obtained  between  measured  and  predicted  acoustic  sig¬ 
natures  using  the  mean  nondimens ional  force  constants  for  vortex  shed¬ 
ding  extracted  from  the  data  analysis,  it  is  believed  that  the  basis  of 
a  realistic  method  of  predicting  the  total  acoustic  signature  of  any 
helicopter  rotor  in  various  flight  regimes  has  been  developed  and 
demonstrated  . _ _ _ _  _ 


l|fl|  4  J  ■« A  M^LACIB  OD  POIIH  1471.  I  jAN  14,  WHICH  ll 

*  m7m)4  i  A  r°"  amv  uM 


Unclassified 

Security  ClMilflcatlofl 


DEPARTMENT  OF  THE  ARMY 

U.  S.  ARMY  AIR  MOBILITY  RESEARCH  *  DEVELOPMENT  LABORATORY 

IUSTIS  DIRECTORATE 
PORT  EU8TI8.  VIRGINIA  23804 


This  report  has  been  reviewed  by  the  Eustls 
Directorate,  U.  S.  Army  Air  Mobility  Research 
and  Development  Laboratory  and  is  considered 
to  be  technically  sound. 

This  program  was  initiated  to  study  helicopter 
rotor  vortex  noise  both  analytically  and  ex¬ 
perimentally.  In  pursuing  this  objective,  it 
was  necessary  to  make  simplifying  assumptions 
to  produce  a  tractable  problem.  Some  of  these 
assumptions,  particularly  those  related  to 
Strouhal  number,  are  acknowledged  to  represent 
limitations  in  the  analysis.  However,  these 
limitations  are  dictated  by  the  state  of  the 
art  and  should  be  treated  as  fertile  areas  for 
future  research,  not  as  errors. 

Generally,  this  report  documents  the  results 
of  unique  and  pertinent  research  relative  to 
a  source  of  rotor  noise  about  which  little 
was  known  prior  to  this  work.  The  report  is 
published  for  dissemination  of  this  data  and 
the  stimulation  of  ideas  relative  to  the  subject. 

The  program  was  conducted  under  the  technical 
management  of  CPT  T.  D.  Evans  and  Mr.  W.  E. 
Nettles  of  the  Aeromechanics  Division  of  this 
Directorate. 


Task  1F162204AA4104 
Contract  DAAJ02-70-C-0023 
USAAMRDL  Technical  Report  72-2 
February  1972 


INVESTIGATION  OF  THE  VORTEX  NOISE 
PRODUCED  BY  A  HELICOPTER  ROTOR 

Final  Report 

RASA  Report  71-10 


By 

H.  Kevin  Johnson  and  Walter  M.  Katz 


Prepared  by 

Rochester  Applied  Science  Associates,  Inc. 
Rochester,  New  York 


for 

EUSTIS  DIRECTORATE 

.  S.  ARMY  AIR  MOBILITY  RESEARCH  &  DEVELOPMENT 

LABORATORY 

FORT  EUSTIS,  VIRGINIA 


Approved  for  public  release; 
distribution  unlimited. 


ABSTRACT 


Karman-street-type  vortex  shedding  from  a  lifting  surface  was 
analyzed  as  a  source  of  noise  from  a  helicopter  rotor  in  hover 
and  forward  flight.  The  theoretical  program  that  was  conducted 
used  acoustic  pressure-time  histories  measured  by  NASA  at  the 
Wallops  Island  Test  Station  for  a  UH-1B  helicopter  in  hover  and 
forward  flight. 

The  experimental  pressure-time  histories  were  analyzed  and  high 
resolution  spectra  were  developed  over  a  frequency  range  of  0  to 
5000  Hz  using  a  0.7-Hz  filter.  On  the  basis  of  these  spectra, 
the  main  and  tail  rotor  rotational  noise  and  discrete  noise 
sources  were  identified  and  then  removed  from  the  measured 
pressure-time  histories,  leaving  pressure-time  histories  repre¬ 
senting  only  the  broadband  noise  radiated  from  the  rotor  system. 
This  broadband  noise  was  then  semiempirically  related  to  the 
noise  induced  by  Karman-vortex  street  shedding  on  the  rotor 
blade. 

The  theoretical  analysis  that  was  developed  related  the  acoustic 
radiation  of  the  vortex  shedding  forces  on  the  blade  to  the 
pressure-time  histories  measured  by  the  microphones  located  in 
the  far  field.  This  analysis  was  then  used  to  solve  for  the 
oscillatory  forces  on  the  blade  that  would  duplicate  the  broad¬ 
band  noise  characteristics  of  measured  pressure-time  histories 
of  the  helicopter  in  various  flight  regimes. 

The  results  of  the  investigation  indicated  that  "vortex  noise" 
is  the  major  source  of  acoustic  radiation  from  a  helicopter 
rotor  in  hover  or  low-speed  flight  and  that  it  is  concentrated 
in  the  frequency  range  of  200  to  500  Hz.  Because  of  the 
excellent  correlation  obtained  between  measured  and  predicted 
acoustic  signatures  using  the  averaged  nondimens ional  force 
constants  for  vortex  shedding  extracted  from  the  data  analysis, 
it  is  believed  that  the  basis  of  a  realistic  method  of  predict¬ 
ing  the  total  acoustic  signature  of  any  helicopter  rotor  in 
various  flight  regimes  has  been  developed  and  demonstrated. 
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This  program  was  conducted  by  Rochester  Applied  Science 
Associates,  Inc.,  under  Contract  DAAJ02-70-C-0023, 

Task  1F162204AA4104 ,  and  was  carried  out  under  the 
technical  cognizance  of  CPT  Timothy  D.  Evans  and 
Mr.  William  E.  Nettles,  Eustis  Directorate,  U.  S.  Army 
Air  Mobility  Research  and  Development  Laboratory, 

Fort  Eustis,  Virginia. 

The  principal  investigator  at  RASA  was  Dr.  H.  Kevin  Johnson; 
Dr.  Walter  M.  Katz  was  responsible  for  the  numerical  analysis; 
Mr.  Lawrence  R.  Sutton  was  involved  with  several  phases  of 
the  program.  Boeing-Vertol  generated  the  spectra  with  a 
Ubiquitous  Spectrum  Analyzer. 
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INTRODUCTION 


The  primary  contributor  to  the  external  sound  of  modern  gas- 
turbine  powered  helicopters  is  the  aerodynamically  induced 
noise  from  the  rotor  system.  Noise  measurements  have  shown 
that  the  acoustic  pressure-time  history  at  an  observer's  loca¬ 
tion  is  due  almost  entirely  to  the  noise  output  of  main  and 
tail  rotor  systems.  This  means  that  the  helicopter's  detect¬ 
ability,  and  to  a  large  extent  its  effectiveness  as  a  weapons 
system,  is  determined  by  the  noise  signature  of  its  rotor 
system.  In  addition,  because  of  the  strong  emphasis  today  in 
noise  pollution  control,  there  has  been  interest  in  rotor  noise 
with  regard  to  commercial  helicopter  operations.  For  these 
reasons,  aerodynamically  induced  sound  from  rotor  systems  has 
been  extensively  investigated,  particularly  in  the  last  5  years. 
(See  References  1,  2,  3,  4,  and  5.) 

Rotor  noise  consists  of  two  basic  types  of  sound  signatures:  one 
is  repetitive  in  frequency;  the  other  is  broadband  or  nonrepeti- 
tive.  The  repetitive  noise  is  typically  called  "rotational" 
noise  in  which  the  frequencies  are  integral  multiples  of  the 
rotor-blade  passage  frequency.  The  broadband  or  ncnrepetitive 
noise  can  be  generally  classified  as  "vortex"  noise.  In  addi¬ 
tion,  there  are  other  helicopter  noise  classifications  such  as 
"blade  slap",  and  sources  such  as  engine  and  transmission  noise. 
These  other  noises  and  sources  are  not  treated  explicitly  in 
this  report.  (See  Reference  5  for  a  discussion  of  rotor  noise 
sources . ) 

A  standard  method  by  which  to  reduce  acoustic  data  for  analysis 
is  to  generate  plots  of  the  noise  spectra.  These  represent  the 
acoustic  energy  distribution  as  a  function  of  frequency.  It  is 
well-known  that  spectra  generated  from  the  same  data  by  differ¬ 
ent  sound-analyzing  devices  can  vary  significantly.  This  is 
partly  due  to  the  fact  that  the  spectra  are  dependent  upon  the 
length  of  record  sampled  and  the  effective  filter  bandwidth  of 
the  device.  In  addition,  the  length  of  record  selected  estab¬ 
lishes  a  lower  bound  on  the  bandwidth  resolution  that  can  be 
obtained.  This  bound  is  independent  of  the  bandwidth  of  the 
analyzing  device. 

An  example  of  the  variation  in  spectra  generated  from  the  same 
data  is  shown  in  Figures  1  and  2.  Shown  are  spectra  developed 
using  a  Ubiquitous  Analyzer  (bandwidth  *  2  Hz;  length  of 
record  =  0.5  sec)  and  a  Bruel  &  Kjaer  Analyzer  with  constant 
percentage  bandwidth  (bandwidth  =  6%) .  One  section  of  this 
report  will  explore  these  differences.  It  is  believed  that  a 
study  of  the  effects  of  filter  bandwidth  and  the  length  of 
record  on  the  noise  spectra  is  important  for  proper  analysis  of 
helicopter  noise  by  this  means. 
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Rotational  noise  is  typically  regarded  as  dominant  for  rotor 
systems  (see  Lowson,  Reference  5).  For  this  reason,  many 
investigators  have  concentrated  in  this  area  of  acoustic  re¬ 
search.  Results  have  shown  that  whereas  rotational  noise  is 
dominant  in  most  cases  for  frequencies  below  100  Hz,  above  100 
Hz,  rotational  noise  may  dominate,  but  in  most  cases  broadband 
noise  will  mask  out  the  rotational  noise.  It  has  been  found, 
however,  that  the  dominant  noise  is  dependent  upon  the  location 
of  the  noise  source  relative  to  the  observer.  For  example, 
refer  to  Figures  1  and  3.  These  spectra  were  created  from  the 
same  data  source  and  illustrate  the  change  in  contribution  of 
the  various  sources  of  noise  to  the  acoustic  signature  with 
change  in  observer  orientation  and  distance.  The  noise  data 
were  taken  for  a  UH-1B  helicopter  hovering  at  100  feet  altitude. 
Figure  3  shows  that  rotational  noise  is  dominant.  This  is  the 
spectrum  developed  from  data  received  at  a  position  700  feet  to 
the  right  side  of  the  helicopter.  In  contrast,  Figure  1  shows 
broadband  noise  to  be  dominant.  This  data  was  taken  at  a  posi¬ 
tion  200  feet  to  the  right  of  the  same  helicopter. 

Since  in  general,  broadband  noise  will  be  nonperiodic  with  re¬ 
spect  to  blade  passage  frequency,  such  noise  is  often  classified 
as  nonrotational .  The  broadband  noise  that  is  associated  with 
helicopter  rotor  blades  can  be  related  to  vortex  shedding  that 
occurs  along  the  blade.  For  the  purposes  of  this  study,  this 
definition  will  be  accepted  and  used  in  the  ensuing  analysis. 
Since  vortex-generated  noise  is  believed  to  occur  in  the  fre¬ 
quency  range  closer  to  that  of  maximum  human  ear  sensitivity, 
the  subjective  loudness  of  such  noise  can  be  greater  than  that 
of  low  frequency  rotor  rotational  noise,  even  though  the  latter 
may  be  of  higher  absolute  magnitude. 

Little  is  known  about  vortex  street  shedding  from  lifting  sur¬ 
faces.  Vortex  shedding  noise  from  rotating  cylindrical  bodies 
was  first  measured  by  Stowell  and  Deming  (6) .  Later  Yudin  (7) 
and  Blokhintsev  (8)  developed  theories  predicting  such  noise. 
Hubb?rd  ( 9  )  has  reported  overall  noise  levels  due  to  all 
sour< es  on  propeller  blades,  and  has  given  an  estimate  of 
noise  due  to  vortex  shedding.  He  found  for  propellers  that 
the  rotational  and  vortex  components  can  be  of  the  same  oru  .. 
of  magnitude.  He  noted  that  the  vortex  component  has  a  higher 
frequency  content  than  the  rotational  component  and  increases 
in  intensity  with  increasing  tip  speed.  Schlegel,  et  al  (3) 
have  made  further  refinements  in  Hubbard's  vortex  noise  magni¬ 
tude  formula.  Krzywoblocki  (10)  has  measured  the  vortex  shed¬ 
ding  frequency  for  airfoils  at  low  Reynolds  numbers;  however, 
the  data  for  angles  of  attack  below  stall  had  wide  scatter. 
Roshko  (ll)  measured  the  frequency  of  shedding  from  circular 
cylinders  in  high  Reynolds  number  flow,  but  not  for  airfoil- 
type  sections. 
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To  summarize  the  experimental  work  to  date,  measurements  have 
been  made  of  the  frequency  of  vortex  shedding,  but  not  within 
the  range  of  flow  conditions  encountered  by  helicopter  rotor 
blades.  To  date,  no  measurements  of  the  oscillatory  lift  and 
drag  forces  due  to  vortex  shedding  have  been  made  within  any 
range  of  flow  conditions.  It  is  therefore  necessary  that  the 
basic  characteristics  of  vortex  shedding  from  airfoils  be  docu¬ 
mented  in  order  to  develop  a  better  understanding  of  their 
effects  on  rotor  noise.  Empirical  or  semiempirical  constants 
will  therefore  have  to  be  used  to  investigate  the  relative 
importance  of  vortex  shedding  on  the  noise  characteristics 
helicopter  rotors  until  suitable  measurements  are  obtained. 

Sadler  and  Loewy  (2)  have  determined,  by  semiempirical  means, 
overall  constants  for  the  lift  and  drag  forces  associated  with 
discrete  vortex  shedding.  In  this  case  the  constants  were 
assumed  to  be  independent  of  Mach  number,  angle  of  attack, 
blade  radius,  and  azimuth.  However,  these  force  constants  may 
not  be  entirely  independent  of  the  aforementioned  parameters. 

The  program  reported  on  herein  was  set  up  so  that  the  lift  and 
drag  force  constants  might  be  developed  from  experimental  data. 
For  the  mathematical  model,  the  helicopter  rotor  blade  was  repre¬ 
sented  acoustically  by  a  series  of  dipoles  whose  strength  and 
frequency  varied  radially  in  the  case  of  hover,  or  both  radi¬ 
ally  and  azimuthally  in  the  case  of  forward  flight.  If  the 
local  angle  of  attack  and  Mach  number  are  known  at  a  blade  ele¬ 
ment,  then  based  upon  Strouhal  number,  the  frequency  of  vortex 
shedding  could  be  precalculated.  Since  there  is  a  known  rela¬ 
tionship  between  dipole  strength  at  a  given  point  in  space  and 
acoustic  pressure  at  another  point,  it  was  possible  in  this 
program  to  use  experimental  pressure-time  histories  measured  at 
a  fixed  observer's  location  to  determine  the  vortex  constants 
that  produced  the  nonharmonic  (broadband)  characteristics  of 
the  measured  signature. 

Therefore,  the  twofold  purpose  of  the  program  discussed  herein 
was  to:  (1)  use  measured  noise  pr^  ■'sure-time  histories  in  con¬ 
junction  with  theoretical  analyses  to  determine  the  effect  of 
bandpass  filter  characteristics  of  plots  of  sound  pressure  level 
versus  frequency;  and  (2)  determine  the  characteristics  of  the 
vortex  shedding  forces  on  a  helicopter  rotor  blade  in  hover  and 
forward  flight. 
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ANALYSIS  OF  EXPERIMENTAL  DATA 


TEST  PROGRAM 

The  data  used  in  the  analysis  were  generated  during  controlled 
acoustic  tests  of  a  UH-1B  helicopter  conducted  during  November 
1969  at  Wallops  Island,  Virginia  (12).  While  these  data  have 
been  discussed  by  Evans  and  Nettles  (12) ,  the  pertinent  aspects 
of  these  tests  will  be  reviewed  here. 

Prior  to  development  of  this  program,  suitably  well-controlled 
acoustic  data  for  helicopters  did  not  exist.  It  was  therefore 
recommended  that  the  Army  conduct  these  tests  with  special  refer 
ence  to  certain  requirements  regarding  documentation  of  test  con 
ditions  and  flight  parameters. 

The  noise  tests  were  conducted  by  the  Army  at  the  Wallops  Island 
Air  Station  with  cooperation  of  the  Dynamic  Loads  Division  of 
the  NASA/Langley  Research  Center.  The  test  area  consisted  of 
an  array  of  microphones  in  the  form  of  a  cross.  This  is  shown 
in  Figure  4.  As  shown,  the  positive  X-axis  is  in  the  direction 
of  flight  and  microphones  1,  2,  and  3  are  at  X  =  200,  500,  and 
700  feet;  microphones  7,  8,  and  9  are  at  X  =  -  200,  -  500,  and 
-  700  feet;  microphones  4,  5,  and  6  are  to  the  right  of  the 
flight  path  at  Y  =  -  200,  -  500,  and  -  700  feet;  microphones  10, 
11,  and  12  are  to  the  left  of  the  flight  path  at  Y  =  200,  500, 
and  700  feet;  and  microphone  13  is  at  X  =  Y  =  0.  A  local  coor¬ 
dinate  system  centered  at  the  main  rotor  hub  moves  with  the 
helicopter  and  is  shown  in  Figure  5.  A  variety  of  test  condi¬ 
tions  were  covered  during  the  test  program,  including  low, 
medium  and  high-speed  forward  flight  and  two  hover  flight  con¬ 
ditions  (6  feet  and  100  feet);  see  Table  I. 

Measured  noise  signatures  were  recorded  on  a  seven-channel 
1/2-inch  tape  along  with  proper  calibration  information.  The 
system  response  was  flat  from  20  to  5000  Hz.  Noise  signatures 
were  also  graphically  recorded  on  an  oscillograph.  The  heli¬ 
copter  position  was  determined  using  radar  and  grapnical  dis¬ 
play  with  time  codes  in  order  to  relate  to  the  acoustic  signal. 
However,  this  means  of  positioning  did  not  meet  the  requirements 
of  the  analysis  to  determine  the  characteristics  of  vortex 
shedding,  as  will  be  discussed  later.  Instead,  the  measured 
acoustic  data  by  itself  could  be  used  to  obtain  the  more  accu¬ 
rate  positioning  of  the  aircraft  as  required  by  the  analysis 
procedures . 

Other  shortcomings  of  the  test  data  for  purposes  of  the  program 
are  as  follows:  helicopter  rotor  rotational  speed  was  not  moni¬ 
tored  accurately,  which  meant  that  small  shifts  in  rotor  rpm  had 
to  be  estimated  from  analysis  of  the  data.  Fuel  consumption  had 
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TABLE  I.  FLIGHT  CONDITIONS  FOR  WHICH  NOISE  SIGNATURES 
WERE  RECORDED  AT  WALLOPS  ISLAND 


Test 

Number 

Flight 
Velocity (kn) 

Altitude (ft) 

Remarks 

1 

54 

100 

- 

2 

54 

100 

- 

3 

83 

100 

- 

4 

83 

100 

- 

5 

115 

100 

Lost  Radar 
Lock -On 

6 

115 

100 

- 

7 

115 

100 

- 

8 

0 

100 

OGE  Hover 

9 

0 

100 

OGE  Hover 

10 

0 

6 

IGE  Hover 

11 

0 

6 

IGE  Hover 

12 

0 

6 

IGE  Hover 

13 


0 


6 


IGE  Hover 


not  been  monitored,  so  it  was  not  possible  to  determine  the  exact 
aircraft  weight  at  a  given  instant  of  time.  Since  rotor  thrust, 
which  in  turn  was  related  to  blade  element  angle  of  attack,  could 
not  be  precisely  determined,  these  values  instead  had  to  be 
estimated  for  the  analysis. 

Perhaps  the  most  desired  feature  not  provided  in  the  data  was 
that  the  UH-1B  rotor  system  had  not  been  instrumented  so  that 
blade  azimuthal  position  could  be  determined  relative  to  the 
recorded  acoustic  signal.  This  limitation  in  the  data  was  over¬ 
come  during  the  analysis  phase  of  the  program  by  making  use  of 
the  measured  rotational  noise  signal.  The  rotor  system  was  set 
at  several  assumed  positions  in  retarded  time,  and  the  rota¬ 
tional  noise  signature  at  an  observer's  location  was  calculated 
for  each  of  these  positions.  By  comparing  phasing  of  the  com¬ 
puted  noise  signatures  with  that  measured,  it  was  then  possible 
to  estimate  true  rotor  blade  position  in  retarded  time. 

DATA  ANALYSTS 


General  Data  Characteristics 


Selected  oscillograph  records  from  the  data  taken  at  Wallops 
Island  are  shown  in  Figures  6  through  8.  Presented  are  time 
histories  of  the  pressure  recorded  at  a  given  microphone  loca¬ 
tion.  The  scales  shown  were  established  from  the  recorded 
calibration  signal.  A  nontrivial  problem  in  reducing  the 
acoustic  data  is  that  of  determining  polarity  or  the  direction 
of  positive  pressure.  This  was  found  by  analytical  means  by 
relating  a  moving  dipole  to  the  measured  pressure-time  histo¬ 
ries.  Since  the  figures  in  this  report  are  plotted  as  they 
were  recorded,  positive  pressure  (where  positive  denotes  a  com¬ 
pressive  sense)  is  the  reverse  of  that  indicated  on  the  plots. 

Figure  6  shows  the  pressure-time  history  recorded  at  position 
4  with  the  helicopter  in  the  100-ft  hover  configuration  over 
the  data  cross.  The  large  peaks  that  are  sinusoidal  in  nature 
occur  at  the  blade  passage  frequency  and  are  associated  with 
rotational  noise.  The  high  frequency  noise  superimposed  on  the 
rotational  noise  is  assumed  to  be  broadband  or  "vortex”  noise. 
This  noise  is  not  periodic  with  blade  passage  frequency.  The 
tail  rotor  rotational  noise  signature  is  largely  suppressed  in 
this  trace  because  the  plane  of  the  tail  rotor  is  nearly  per¬ 
pendicular  to  a  line  drawn  from  the  recording  microphone  at 
position  4  to  the  tail  rotor. 

As  expected,  the  forward  flight  data  differ  from  those  taken 
for  hover,  even  at  the  same  microphone  position.  The  effect 
of  distance  and  directivity  is  shown  by  comparing  Figure  7, 
taken  early  in  the  flyby  record,  with  Figure  8,  taken  when  the 
helicopter  was  near  ground  zero  of  the  data  cross.  In  the  first 
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figure,  where  the  helicopter  is  at  a  distance,  the  tail  rotor 
noise  peaks  are  clearly  evident  and  occur  at  about  five  times 
the  frequency  of  the  main  rotor  peak.  (The  ratio  of  tail-to- 
main  rotor  rotational  speed  is  5.108  to  1.)  In  Figure  8,  where 
the  helicopter  is  near  the  observer,  tail  rotor  rotational  noise 
does  not  show  as  clearly.  The  difference  is  due  to  the  change 
in  directivity  in  that  when  the  helicopter  is  far  away  and 
approaching  the  observer,  the  observer  is  essentially  in  the 
plane  of  both  rotors.  Conversely,  when  the  helicopter  is  near 
the  microphones,  the  plane  of  the  tail  rotor  is  approximately 
perpendicular  to  a  line  drawn  from  position  4  to  the  tail  rotor. 
Also,  the  broadband  or  vortex  noise  is  not  evident  when  the  heli¬ 
copter  is  very  far  away  but  becomes  increasingly  important  as  the 
helicopter  approaches  the  microphones.  (See  Figures  7  and  8.) 
This  effect  is  also  apparent  in  the  spectra  shown  in  Figures  9, 
10,  and  11,  which  represent  spectra  for  one  flyby  taken  respec¬ 
tively  with  helicopter  position  3000  ft  from  ground  2ero 
(approaching),  at  ground  zero,  and  at  1500  ft  from  ground  zero 
(moving  away). 

Doppler  Effect 

The  Doppler  effect  can  be  seen  from  the  forward  flight  oscillo¬ 
graph  traces.  The  change  in  frequency  that  occurs  during  the 
flyby  is  observed  by  comparing  the  distances  on  the  plots  (times) 
between  tail  rotor  peaks.  See  Figures  12  and  13.  The  shorter 
distances  or  times  are  shown  in  Figure  12  and  correspond  to  a 
higher  frequency  than  those  shown  in  Figure  13.  The  Doppler 
shift  depends  on  the  helicopter's  velocity,  direction  of  flight 
and  its  instantaneous  position  with  respect  to  the  observer. 

For  example,  assume  that  an  observer  is  located  at  ground  zero 
and  the  helicopter  is  approaching  the  observer  head-on.  The 
apparent  frequency  of  a  pure-tone  sound  source  on  the  helicopter 
would  be  increased  by  a  constant  amount  during  the  approach; 
then  there  would  be  a  steep  drop  in  frequency  as  the  helicopter 
passed  ground  zero.  As  the  hel.  copter  moved  away  from  the 
observer,  the  apparent  frequency  would  be  decreased  by  a  con¬ 
stant  amount.  For  the  measured  ita,  the  observer  (microphone) 
is  off  to  one  side  of  the  flight  >ath.  Decause  of  this,  when 
the  helicopter  is  in  the  distance  the  time  interval  between 
tail  rotor  peaks  remains  relative^  -  constant,  corresponding 
somewhat  to  the  helicopter  approac  ing  the  observer  head-on 
case.  On  the  other  hand,  as  the  ht  icopter  gets  nearer  the 
observer,  there  is  now  a  gradual  dro,  in  frequency  from  the 
initial  increased  value,  rather  than  >.  e  steep  drop  as  xn  the 
head-on  case.  This  was  evident  from  ti.  records  where  it  war. 
noted  that  as  the  helicopter  approached  *  ->arer,  the  time  between 
tail  rotor  peaks  increased,  corresponding  to  a  decrease  in  fre¬ 
quency  with  respect  to  the  initial  increased  value. 
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The  sound  from  the  blade  in  an  advancing  position  relative  to 
an  observer  dominates  the  noise  generated  from  the  blade  in  all 
other  azimuth  positions.  This  fact  makes  it  easier  to  under¬ 
stand  certain  aspects  of  the  data  that  were  used  herein.  For 
example,  the  pressures  recorded  from  hover  at  microphones  6  and 
12,  700  ft  to  the  right  and  left  of  the  flight  path  respec¬ 
tively,  are  not  the  same.  See  Figures  14  and  15.  This  differ¬ 
ence  is  due  in  part  to  asymmetry  associated  with  the  helicopter. 
For  example,  in  Figure  14,  taken  to  the  right  of  the  helicopter, 
a  portion  of  the  tail  rotor  disk  is  shielded  from  observers  to 
the  right  by  the  tail  pylon.  Observers  to  the  left  (Figure  15) 
are  not  so  shielded.  The  net  effect  of  this  shielding  is  to 
make  the  pressure-time  histories  differ  as  received  by  observers 
on  eithor  side  of  the  flight  path.  Other  asymmetries  noted  in 
the  data  were  caused  by  the  helicopter  not  holding  position  and 
orientation.  Since  no  time  code  was  recorded  during  the  hovers, 
the  record  from  one  microphone  position  could  not  be  correlated 
to  that  of  another  microphone  position.  Also,  the  advancing 
blade  as  seen  from  position  6  passes  near  the  wake  of  the  tail 
rotor  and  across  the  tail  boom,  while  the  advancing  blade  as 
seen  from  position  12  passes  over  the  nose  of  the  helicopter 
through  relatively  undisturbed  air.  Tnis  effect  is  significant 
because  the  sound  of  the  advancing  blade  dominates  the  noise 
from  the  blade  in  all  other  positions  relative  to  an  observer. 
Finally,  asymmetries  caused  by  the  wind  blowing  across  tho  flight 
path  in  varying  direction  and  magnitude  can  distort  tho  sound 
signature.  For  these  tests,  wind  speed  varied  from  2  to  9  knots 
with  a  mean  velocity  of  about  5  knots,  and  the  direction  changed 
from  cast  to  south  with  a  mean  of  about  150*.  In  the  case  of 
tho  hover  data,  there  were  no  records  available  to  indicate  tho 
precise  position  of  tho  helicopter. 

The  change  in  tho  character  of  the  tail  rotor  noise  with  time 
during  an  approach  is  shown  in  Figure  16.  Observe  tho  splitting 
of  the  tail  rotor  pulses  as  the  aircraft  approaches  ground  zero. 
This  splitting  of  the  pulse  has  been  a  subject  of  controversy, 
but  can  bo  properly  explained  as  a  ground  reflection  effect. 

The  earlier  time  that  is  shown  in  the  record  corresponds  to  the 
helicopter  being  further  away.  For  a  fixed  microphone  height, 
increasing  the  distance  between  a  noise  source  and  the  observer 
will  decrease  the  time  delay  between  the  original  and  ground 
reflected  signal.  As  the  helicopter  approaches  the  observer, 
the  reflected  signal  of  the  tail  rotor  becomes  more  and  more 
retarded  from  the  original  signal,  until  two  distinct  signals 
result.  This  means  that  the  earlier  pulse  signals  shown  in 
thesQ  records  actually  represent  a  superposition  of  the 
original  and  reflected  signals. 
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Ground  Reflection  Effect 


The  importance  of  ground  reflection  in  the  data  is  revealed  by 
comparing  spectra  for  different  microphone-helicopter  orienta¬ 
tions;  see  Figures  17  and  18,  which  represent  extended  spectra 
at  positions  4  and  6.  The  general  characteristic  of  each  spec¬ 
trum  $  n  the  broadband  region  reveals  a  series  of  maxima  and 
minirv  separated  throughout  the  record  by  equal  frequency  inter¬ 
vals.  However,  this  frequency  interval  is  not  the  same  for 
both  recording  positions.  The  frequency  interval  between  suc¬ 
cessive  maxima  or  minima  on  the  records  is  seen  to  change  as  the 
orientation  between  source  and  observer  changes.  For  a  given 
microphone  height  above  the  ground,  the  microphone  receives  one 
signal  directly  from  the  sound  source  and  another  that  is  re¬ 
flected  from  the  ground  and  then  received.  The  frequency  in¬ 
terval  referred  to  between  successive  peaks  or  valleys  on  the 
SPL  records  is  related  to  this  time  separation  between  original 
and  reflected  signal.  Thus,  the  frequency  spectrum  received  by 
an  obsorver  is  distorted  from  that  of  the  source;  some  frequen¬ 
cies  will  be  amplified  by  reflection  while  others  arc  suppressed. 

As  noted,  the  acoustic  signature  at  a  microphone  (from  a  point 
source)  i3  a  linear  combination  of  the  directly  incident  sound 
wave  and  the  reflected  sound  wave.  Assuming  specular  reflec¬ 
tion,  the  reflected  sound  wave  i3  proportional  to  the  initial 
acoustic  signature  delayed  by  the  time  needed  for  the  sound 
wavs  to  propagate  the  longer  distance. 

For  a  source  of  altitude  h,  a  distance  i  from  a  microphone 
located  a  height  «  off  the  ground,  the  additional  distance  A 
traveled  by  the  reflected  wave  is 


a  -  /(b+r )  -  /Vo  2+&2  (1) 

When  t  is  small  compared  to  {j.  and  A  is  well  approximated  by 
the  expression 


(2) 


The  time  delay  the  reflected  wave  is  therefore 

T  - 

c*42+d 


(3) 


whore  c  is  the  speed  of  sound. 
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For  a  given  distance  and  microphone  height,  when  the  frequency 
at  the  microphone  of  the  reflected  sound  wave  is  shifted  by 
180°  or  multiple  thereof  from  that  of  the  direct  sound  wave, 
the  sound  waves  will  interfere  and  cancel.  This  accounts  for 
the  minima  in  the  SPL  plots  shown.  In  contrast,  if  the  phase 
shift  between  direct  and  reflected  waves  is  360*  or  multiple 
thereof,  then  the  incoming  sound  signal  at  t <£*  TnliVophone  will 
be  reinforced  due  to  direct  addition  of  ^he  tw  sisals.  This 
is  the  explanation  for  the  maxima  on  the  SPa, 

The  frequencies  which  correspond  to  a  iBO*  a#  shift  are 


f  *  c  (n  ♦  1/2)  (4) 

where  n  is  a  non-negative  integer.  These  f rt«.jliencie*  cancel. 
For  the  frequencies 


the  direct  and  reflected  waves  are  in  phase  and  add. 

This  effect  is  dramatically  demonstrated  in  the  spectra  of  the 
helicopter  in  a  100  ft  hover,  as  shown  in  Figures  17  and  18  for 
positions  4  and  6,  respectively.  (These  positions  arc  200  and 
700  feet  to  the  side,  both  microphones  are  5.5  feet  high.) 

The  reflected  signal  time  lag  at  position  6  is  less  t-nan  that 
of  position  4  because  the  increase  in  distance  trav  led  by  a 
reflected  wave  received  at  position  6  is  less  than  .  a  increase 
in  distanco  travelled  by  a  reflected  wave  received  at  position 
4.  That  is,  the  increase  in  distance  travelled  by  a  reflected 
wave  at  position  6  is  1.48  ft,  while  the  increase  in  distance 
travelled  by  a  reflected  wave  at  position  4  is  4.65  ft.  There¬ 
fore,  the  time  lag  for  a  sound  wave  at  position  6  is  1.34  x 

10”  sec  (745  Hz),  while  the  tine  lag  for  position  4  is  4.22  x 

10”  sec  (237  Hz).  The  maxima  or  reinforced  frequencies  would 
therefore  correspond  to  integral  multiples  of  this  frequency, 
and  the  minima  or  cancelled  frequencies  would  correspond  to 
(n  ♦  1/2) f.  Tho  first  minimum  frequency  at  position  6  is  373 
Hz,  and  the  first  and  second  minima  at  position  4  arc  118  liz 
and  356  Hz,  respectively  (see  Figures  17,  18).  Thoso  frequen¬ 
cies  arc  consistent  with  the  predicted  minima  (372.5  Hz  for 
position  6  and  118.5  Hz  and  355.5  Hz  for  position  4). 

The  first  minima  of  positions  4  and  10,  which  are  equidistant 
from  the  data  cross,  did  not  occur  at  the  same  frequency.  This 
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indicated  that  the  helicopter  was  not  in  precise  position.  In 
fact,  calculations  showed  that  the  helicopter  was  only  179  ft 
from  position  4  and  95  ft  high,  not  at  200  ft  and  100  ft  high. 
With  these  distances  substituted  into  the  previous  calculations, 
the  calculated  maxima  and  minima  checked  closely. 

Digitization 

In  order  to  perform  computer  analysis  of  the  acoustic  data,  the 
pressure-time  history  had  to  be  converted  into  digital  form  to 
be  compatible  with  computer  input  format.  This  process,  called 
digitization,  was  first  performed  manually  on  oscillograph 
records  which  had  been  made  directly  from  the  original  FM  record¬ 
ings.  The  hover  oscillograph  trace  was  digitized  graphically 
with  an  800  Hz  calibration  signal  as  a  reference.  The  hover 
record  digitized  in  this  manner  was  used  for  the  filter  bandwidth 
analysis  phase  of  the  contract.  Figure  19  compares  the  manually 
digitized  data  with  the  oscillograph  record. 

However,  while  the  frequency  content  of  the  manually  digitized 
data  was  adequate  for  the  filter  analysis,  it  was  inadequate  for 
carrying  out  the  vortex  shedding  analysis.  The  original  FM 
recordings  were  then  redigitized  using  a  high-resolution  elec¬ 
tronic  conversion  technique  specially  developed  by  RASA  for 
carrying  out  this  work.  By  coupling  a  computer  with  special 
circuitry,  a  1.45-second  record  could  be  digitized  at  an  11.3 
KHz  rate  providing  approximately  16000  data  points.  The  result¬ 
ing  frequency  response  of  the  digitized  data  extended  from  about 
0.7  Hz  to  5650  Hz,  which  exceeded  the  range  of  the  recording 
system  of  20  Hz  to  5000  Hz.  The  signal-to-noise  ratio  of  the 
developed  system  was  excellent,  at  better  than  45  db. 

To  demonstrate  the  accuracy  of  this  digitization,  tho  digitized 
data  were  plotted  and  compared  with  the  corresponding  oscillo¬ 
graph  records;  seo  Figure  20.  The  digitized  data  are  in  excel¬ 
lent  agreement  with  the  oscillograph  trace.  In  fact,  if  differ¬ 
ences  do  exist,  the  electronically  digitized  results  should  be 
considered  tho  more  reliablo  for  two  reasons.  First,  tho  con¬ 
version  was  done  diroctly  off  the  FM  carrier  of  the  recorded 
signal.  Thib  technique  bypasses  tho  FM  demodulator  unit,  which 
in  itself  introduces  some  distortion  to  the  output.  In  addition, 
the  electronically  digitized  data  bypasses  the  frequency  response 
of  the  oscillograph  galvanometers,  which  were  flat  to  only  about 
3  KC. 
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FILTER  BANDWIDTH  ANALYSIS 


GENERAL  DISCUSSION 


Spectra  generated  from  the  same  pressure-time  history  by  differ¬ 
ent  spectrum  analy2ers  can  vary  significantly  (see  Figures  1  and 
2) .  These  are  spectra  created  from  the  same  data — the  pressure¬ 
time  history  for  a  UH-1B  helicopter  in  a  100-ft  hover  recorded 
at  a  distance  of  200  feet  Figure  1  shows  the  noise  spectrum 
generated  by  a  Ubiquitous  Analyzer  with  an  effective  filter 
bandwidth  of  2  Hz,  whereas  Figure  2  presents  the  spectrum  re¬ 
sulting  from  a  Bruel  6  Kjaer  Analyzer,  which  had  a  constant 
percentage  bandwidth  filter  of  6%.  For  a  constant  percentage 
bandwidth  filter,  the  bandwidth  of  the  filter  increases  with 
frequency,  which  means  that  the  high  frequency  end  of  the  spec¬ 
trum  is  effectively  smeared  out  due  to  the  large  bandwidth. 

This  characteristic  is  well  illustrated  in  Figure  2  where  no 
discrete  peaks  occur  in  the  high  frequency  range,  in  contrast  to 
Figure  1  where  discrete  peaks  are  observed  in  the  high  frequency 
range. 

In  addition  to  the  effects  of  filter  bandwidth,  other  parameters 
influence  the  resulting  spectrum.  Perhaps  the  most  important  of 
these  is  the  length  of  record  over  which  the  spectrum  is  devel¬ 
oped.  Since  spectra  created  from  digital  information  are  in 
effect  plots  of  the  resulting  Fourier  coefficients,  these  spectra 
will  depend  on  the  time  interval  over  which  the  data  is  sampled. 
In  fact,  for  a  Fourier  analysis  in  which  no  filter  characteris¬ 
tics  are  prescribed,  the  effective  bandwidth  of  the  resulting 
spectrum  is  the  reciprocal  of  the  length  of  record.  This  means 
that  a  1-second  record  would  result  in  a  spectrum  plot  limited 
to  a  1  Hz  filter;  a  10-second  record  would  yield  a  spectrum 
with  effectively  a  0.1  Hz  filter,  etc.  The  reciprocal  value  of 
the  record  length  provides  a  lower  bound  on  the  resolution  that 
can  bo  obtainod  in  the  spectrum.  That  is,  no  matter  how  narrow 
the  filtor  bandwidth,  no  more  information  can  bo  determined  from 
a  given  record  than  for  a  filter  of  bandwidth  equal  to  the  recip¬ 
rocal  of  the  length  of  record.  This  means  that  in  general  the 
resolution  will  be  somewhat  less  than  this  optimum  amount  when 
a  filter  is  used  to  develop  a  spectrum. 

Since  the  electronic  equipment  by  which  the  spectra  are  obtained 
is  not  standard,  there  has  frequently  been  controversy  over 
interpretation  of  SPL's  obtained  from  similar  acoustic  pressure¬ 
time  histories.  In  addition,  there  has  been  considerable  dis¬ 
cussion  as  to  what  filter  bandwidth,  record,  length,  etc., 
should  be  used  to  obtain  fully  descriptive  and  accurate  plots 
of  SPL  versus  frequency  for  a  given  sound-pressure  history. 

To  clarify  these  questions,  an  investigation  was  carried  out 
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to  determine  the  effects  of  these  parameters.  The  results  of 
this  investigation  will  be  discussed  in  the  following  paragraphs. 

ANALYTICAL  REPRESENTATION  OF  BANDPASS  FILTER  CHARACTERISTICS 

The  characteristics  of  a  bandpass  filter  are  usually  repre¬ 
sented  by  attenuation  curves  (in  db)  which  show  the  attenuation 
levels  of  signals  of  given  frequencies  relative  to  the  center 
frequency  of  the  filter.  A  representative  attenuation  curve, 
that  of  the  Federal  Scientific  Ubiquitous  Analyzer  used  for 
this  program,  is  presented  in  Figure  21.  In  order  to  apply 
the  attenuation  effects  of  a  filter  to  a  frequency  spectrum* 
it  is  necessary  to  consider  the  attenuation  levels  in  terms  of 
their  related  percentage  factors,  i  .e.*  percentage  of  signal 
amplitude  transmitted  through  filter.  The  conversion  of  atten¬ 
uation  values  to  percentage  values  can  be  achieved  by  the 
expression 

(attenuation) 

“  11  ■■  11  1 

percentage  factor  (in  decimal  notation)  =10  * 

To  obtain  the  percentage  factors  in  a  functional  form,  the 
attenuation  curve  can  he  approximated  by  a  power  series  repre¬ 
sentation  and  substituted  in  the  previous  expression  to  result 
in 


P 


f 


— 

l  a  (AFB)bn 

/20 

n=o 

(6) 


in  which  Pf  is  the  percentage  factor,  an  are  the  polynomial 

coefficients  that  define  the  filter  curve,  AFB  is  a  function  of 
the  differonco  between  the  center  frequency  and  the  frequency 
of  concern,  and  b  is  an  integer  value  of  1  if  the  filter  curve 
is  not  symmetric  about  the  center  frequency  and  2  if  the  curve 
is  symmetric.  A  curve  results  which  has  values  ranging  from 
near  zero  at  the  outer  analytical  cutoff  limits  of  the  filter 
to  one  at  the  center  frequency  point.  For  a  constant  bandwidth- 
type  filter, the  term  AFB  may  be  defined  as 

AFB  (u)  -  u>c ) /BW  (7) 

in  which  <*>  is  the  frequency  of  interest,  u  is  the  center  fro- 

V 

quency  of  the  filter,  and  BW  is  the  filter's  bandwidth. 
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To  carry  out  the  analysis,  a  sat  of  representative  filter  charac¬ 
teristics  had  to  be  selected.  Since  experimental  data  were 
available  from  the  Ubiquitous  Analyzer  by  which  analytical  re¬ 
sults  could  be  compared,  the  filter  characteristics  of  that 
instrument  (see  Figure  21)  were  selected  for  modelling.  The 
Ubiquitous  filter  is  symmetric,  and  hence  b  ■  2  in  Equation  (6). 
The  following  values  of  aR  provided  a  good  fit  to  the  filter 

characteristics  of  the  Ubiquitous  Analyzer: 


a0 

U 

0.0 

ai 

m 

9.07221 

a2 

m 

4.48306 

a3 

m 

-3.28649 

a«, 

m 

0.818003 

a  5 

m 

-0.0904785 

a  6 

s 

0.0037009 

APPLICATION  OF  FILTER  CHARACTERISTICS  TO  DATA 

To  explore  effects  of  filter  bandwidth  and  record  length,  a 
typical  pressure-time  history  recorded  for  the  IJH-1B  in  hover 
was  digitized.  The  digitization  was  carried  out  prior  to 
RASA's  development  of  an  electronic  digitization  technique,  and 
consequently  was  performed  manually.  Each  second  of  record  was 
subdivided  into  1600  intervals  and  pressure  values  were  scaled 
directly  from  the  oscillograph  records.  Therefore,  the  maxi¬ 
mum  frequency  content  for  the  digitized  records  was  800  Hz. 

Digitized  pressuro-time  histories  were  Fourier  analyzed  for 
each  case  in  the  form 

N 

P(t)  ■  l  (A^  cos  u^t  ♦  Bk  sin  w^t)  (8) 

where  A^  and  B^  are  the  Fourier  coefficients  for  each  frequency 

and  N  is  the  number  of  frequencies  which  can  be  obtained 

from  a  given  record  and  depends  on  the  number  of  digitized 
points  of  pressure  versus  time.  Note  that  2N  input  points 
would  be  required  to  obtain  N  sots  of  A^  and  B^. 

The  bandpass  filter  spectra  were  obtained  by  application  of  the 
desired  bandpass  filter  curve  to  the  Fourier  coefficient  infor¬ 
mation  at  the  desired  frequencies.  For  this  analysis,  the 
Ubiquitous  filter  characteristics  with  different  bandwidths 
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were  applied.  To  avoid  unnecessary  numerical  work  for  small 
contributions  from  frequencies  sufficiently  distant  from  the 
center  frequency,  u  ,  an  upper  and  lower  cutoff  frequency  was 

assumed  such  that  only  frequencies  in  the  range 


-  2.5  BW  +  u  “ 
c 


<d. 


2.5  BW  +  u 

c 


(9) 


were  considered  to  contribute  for  a  constant  bandpass  filter. 

For  each  frequency  point  of  the  calculated  spectrum,  the  result¬ 
ing  pressure  was  obtained  by  summing  contributions  of  all  fre¬ 
quencies  that  passed  through  a  set  of  specified  filter  charac¬ 
teristics.  In  the  case  of  a  constant  bandwidth  filter,  set  at 
the  center  frequency,  u  ,  the  pressure  due  to  the  applicable 

w  ^ 

frequencies,  oj^,  near  the  center  frequency  was  found  from 


P(u>  ,BW) 


where 
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and  is  used  to  specify  summation  of  the  effects  of  all 

which  are  in  range  defined  by  Equation  (9) .  The  SPL  for  fre¬ 
quency,  w  ,  is  obtained  by 


SPL(tdc,BW)=201og1  o  [p  (wc  ,BW)/.  0002  j  =201og  ]  0  [p  (idc,BW)/2j +80db 


(11) 

For  calculating  spectra  for  the  case  of  a  constant-percentage- 
type  bandpass  filter, the  above  expressions  should  be  modified 
by  replacing  the  term  BW  with  BWP,  where  BWP  is  the  bandwidth 
defined  such  that  it  is  a  constant  percentage  of  the  center 
frequency,  u  . 

DESCRIPTION  OF  ANALYSIS  CONDUCTED 

Filter  characteristics  of  the  Ubiquitous  Analyzer  where  modelled 
analytically  with  provision  to  treat  different  constant 
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bandwidth  values  as  well  as  constant-percentage-type  filters. 

This  meant  that  while  a  variety  of  widths  of  filter  would  be 
applied ,  each  would  have  the  general  characteristics  of  the 
Ubiquitous  filter  as  shown  in  Figure  21. 

Variations  in  effective  bandwidth  were  studied  in  the  following 
manner.  A  given  record  was  selected.  Chosen  was  the  hover  con¬ 
dition  with  the  UH-1B  at  100  ft  altitude  with  the  pressure-time 
history  as  recorded  at  position  4,  200  feet  from  the  helicopter. 
Four  variations  in  length  of  record  were  selected:  one  blade 
passage  interval  (.0926  second);  two  blade  passage  intervals 
(.1852  second);  five  blade  passage  intervals  (.463  second);  and 
ten  blade  passage  intervals  (.926  second).  The  selected  pro¬ 
cedure  was  to  take  each  of  the  four  lengths  of  record  and 
develop  analytical  SPL  plots  for  each  record  for  the  following 
conditions:  (1)  unfiltered;  (2)  1  Hz  bandwidth  Ubiquitous-type 

filter;  (3)  2  Hz  bandwidth  Ubiquitous-type  filter;  (4)  5  Hz 
bandwidth  Ubiquitous-type  filter;  (5)  10  Hz  bandwidth  Ubiquitous- 
type  filter;  and  (6)  a  constant-percentage  (1%)  Ubiquitous-type 
filter. 

RESULTS  OF  THE  THEORETICAL  ANALYSIS 

Representative  spectra  of  the  four  record  lengths  analyzed  and 
the  six  simulated  analog  filter  widths  are  shown  in  Figures  22 
through  30.  The  effect  of  varying  the  length  of  record  is 
shown  in  Figures  22,  23,  24,  and  25.  Presented  are  spectra 
for  one,  two,  five,  and  ten  blade  passages  in  which  the  spectra 
were  obtained  from  the  basic  Fourier  analyses  without  applying 
specified  filter  characteristics.  This  means  that  for  each  of 
these  records,  the  effective  filter  bandwidth  may  be  considered 
to  be  the  reciprocal  of  the  record  length.  Hence,  Figure  22 
represents  a  spectrum  developed  with  a  10.8  Hz  filter;  Figure 
23  represents  that  for  a  5.4  Hz  filter;  Figure  24*  a  2.2  Hz 
filter;  and  Figure  25* a  1.1  Hz  filter.  Observe  the  increase  in 
resolution  as  record  length  increases  corresponding  to  a  de¬ 
crease  in  filter  bandwidth.  For  example,  the  ten  blade  passage 
record  shown  in  Figure  25  demonstrates  the  discrete  nature  of 
the  main  and  tail  rotor  rotational  noise.  Here,  the  bandwidth 
is  sufficiently  small  to  separate  the  first  harmonic  of  tail 
rotational  noise  (55  Hz)  with  the  fifth  harmonic  of  main  rota¬ 
tional  noise  (54  Hz).  As  record  length  decreases,  the  bandwidth 
increases  commensurately .  The  spectra  in  Figure  24  with  a 
record  length  of  ^-.5  sec  and  a  bandwidth  of  2.2  Hz  can  no 
longer  resolve  the  1  Hz  separation  between  the  first  tail  rotor 
harmonic  with  the  fifth  rotor  harmonic.  When  the  record  length 
is  ^.1  sec* the  bandwidth  is  10.8  Hz  (see  Figure  22).  In  this 
case  the  main  rotor  rotational  harmonics  are  no  longer  resolved. 
The  values  of  the  spectrum  at  these  harmonics  are  equal  to  the 
peak  values  of  the  discrete  harmonic  peaks  in  the  1-sec  record 
shown  in  Figure  25.  These  magnitudes  remain  accurate  because 
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the  length  of  record  for  the  spectrum  in  Figures  22  was  choaan 
to  be  exactly  on*  blade  passage. 

Introducing  a  simulated  filter  of  different  oandwidtha  into  the 
Fourier  analyaia  of  a  given  length  of  racord  has  the  aane  gen-* 
eral  effect  aa  changing  the  length  of  record.  An  increaae  in 
filter  bandwidth  correeponda  to  a  shortening  of  length  of  record. 
Thia  ie  ahown  in  Figurea  26,  27,  28,  and  29.  Shown  are  the 
apectra  from  a  ten-blade-passage  length  of  record  (.926  aecond) 
that  have  been  developed  with  10  Hz,  5  Hz,  2  Hz  and  1  Hz  band¬ 
width  filtera  reepectively.  Observe  the  similarity  between 
these  spectra  and  those  of  Figurea  22,  22,  24,  and  25  in  which 
record  length  was  varied.  For  example,  compare  the  plot  of  the 
one-blade-passage  spectrum  (Fourier  analysis  without  filter) 
shown  in  Figure  22  to  that  of  the  10  Hz  simulated  filter  results 
of  Figure  26.  The  general  features  of  th^se  plots  are  the  sane. 
The  spectrum  in  Figure  26  appears  smoother  than  the  spectrum  in 
Figure  22  because  the  spectrum  in  Figure  20  had  points  calculated 
at  closer  frequency  intervals.  The  structure  of  the  two  plots 
is  about  the  same.  Similar  comparisons  can  he  made  between 
Figures  22,  24,  and  25  and  Figures  7 n ■  2H  Mid  29  respectively. 

Figure  30  illustrates  the  type  of  spectrum  obtained  with  a 
constant-percentage  (1%)  bandpass  filter.  Here,  the  width  of 
the  filter  ia  a  fixed  percentage  of  the  center  frequency.  In 
other  words,  for  a  100  Hz  signal  the  bandwidth  would  be  1  Hs, 
whereas  for  a  1000  Hs  signal  it  woula  be  10  Hz.  Figure  30  ahowa 
that  excellent  resolution  is  obtained  for  Un?  low  frequency  end 
of  the  spectrum  (below  200  Hz)  and  that  this  resolution  deteri¬ 
orates  in  the  higher  frequency  region.  In  the  neighborhood  of 
100  Hz  Figure  30  exhibits  the  same  quality  of  resolution  ob¬ 
tained  throughout  Figure  29  where  a  i-h'z  filter  wan  applied. 
Similarly,  Figure  30  in  the  neighbor  hood  of  2'jc  az  can  be  com¬ 
pared  with  Figure  28  in  which  a  2-Hz  filter  woe*  applied,  etc. 

Using  the  RASA  electronic  digitization  technique  previously 
discussed,  coupled  with  the  Fast  Fourier  Transform,  spectra  of 
significantly  better  quality  can  be  obtained*  nee  Figure  31. 

The  spectrum  shown  in  Figurs  31  if  from  a  15  second  portion  of 
the  hover  record;  the  bandwidth  is  .07  Hz  and  the  frequency 
range  is  0  to  280  Hz.  The  main  and  tail  rotational  noise  har¬ 
monics  appear  to  be  truly  discrete.  Thu  fifth  harmonic  of 
main  rotor  noise  and  the  first  harmonic  cf  tail  rotor  noise 
are  clearly  separated.  In  this  record,  th*  noise  soectrum 
above  100  Hz  is  dominated  by  broadband  rather  then  discrete 
noise.  This  spectrum  also  shows  a  dircrete  peak  at  72  Hz 
which  could  not  be  identified  with  main  or  tai)  rotor  rota¬ 
tional  noise,  nor  with  known  gear  clash  frequencies  of  the 
UH-1B  as  provided  RASA  by  Bell,  interestingly,  this  peak  waa 
also  notad  in  the  Ubiquitous  and  hand-digJ ti zed  spectra  where 
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it  appr  red,  but  not  as  dramatically  as  shown  by  Figure  31  from 
the  RASA  digitized  record.  An  interesting  feature  provided  by 
the  high  resolution  with  RASA's  digitizing  and  spectral  analysis 
procedure  was  the  appearance  of  asymmetry  between  the  two  main 
rotor  blades  in  the  spectrum.  This  is  revealed  in  Figure  31 
by  the  narrow  peaks  10  db  above  background,  at  the  odd  harmon¬ 
ics  of  main  rotor  rotational  speed  (half  harmonics  of  the  blade 
passage  frequency)  which  are  5.4,  16.2,  27  and  37.8  Hz. 


Effects  of  filter  bandwidth  and  recora  length  were  also  con¬ 
sidered  for  records  taken  of  helicopters  in  forward  flight. 
Here,  the  problem  is  complicated  due  to  continual  Doppler 
shift  of  the  frequencies  and  change  of  position  of  the  heli¬ 
copter  with  time.  For  a  sound  source  of  frequency  f  moving 

s 

toward  an  observer  at  velocity  Vg,  the  apparent  frequency  at  a 
listener  is  given  by 


As  an  indication  of  the  magnitude  of  this  effect  consider  a 
helicopter  approaching  a  listener  at  a  velocity  of  100  fps. 
In  that  case 


(13) 


In  other  words,  a  signal  at  100  Hz  would  be  shifted  to  110  Hz, 
that  at  1000  Hz  shifted  to  1100  Hz,  etc. 

This  constant  shift  is  complicated  when  the  observer  is  off  to 
one  side  of  the  flight  path.  In  this  case,  the  frequency  of 
the  sound  source  as  received  by  an  observer  will  be  continually 
changing  with  time.  In  regard  to  a  spectrum  that  might  be 
generated  over  a  record  length  where  this  occurs,  this  means 
that  rather  than  narrow  discrete  peaks  occurring,  the  peaks 
would  instead  be  smeared  and  broadened.  Depending  upon  the 
amount  of  Doppler  shift,  it  would  become  increasingly  difficult 
or  impossible  to  differentiate  discrete  rotational  and  gear- 
clash  noise  sources. 

Another  problem  in  analyzing  forward  flight  records  is  due  to 
the  change  in  orientation  of  the  helicopter  with  time  with 
respect  to  a  fixed  observer.  It  will  be  shown  in  the  next 
section  that  the  primary  contribution  to  rotational  noise  from 
the  rotor  occurs  over  a  very  narrow  azimuth  interval  of  the 
rotor.  This  narrow  azimuth  interval  is  where  the  relative  Mach 
number  of  the  blade  is  the  greatest  with  respect  to  an  observer. 
If,  as  in  the  case  of  forward  flight,  the  helicopter  is  changing 
orientation  with  respect  to  the  observer,  then  the  azimuth 
position  of  the  blade  with  respect  to  its  maximum  contribution 
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to  the  noise  signature  is  also  changing.  At  each  of  these  posi¬ 
tions,  it  will  have  a  different  angle  of  attack  and  velocity 
distribution. 

In  summary,  forward  flight  spectra  are  seen  to  involve  a  trade¬ 
off.  To  minimize  Doppler  frequency  shift  and  changes  in  the 
helicopter's  position,  as  short  a  record  length  as  possible 
should  be  analyzed.  However,  as  record  length  is  reduced,  the 
effective  filter  bandwidth  increases.  This  means  that  if  too 
short  a  forward  flight  record  is  taken,  discrete  noise  sources 
such  as  the  main  and  tail  rotor  rotational  noise  peaks  cannot 
be  distinguished.  On  the  other  hand,  if  record  length  is  in¬ 
creased,  Doppler  shift  and  position  effects  may  similarly  make 
it  difficult  to  distinguish  discrete  frequency  sources.  These 
problem  areas  noted  for  forward  flight  become  more  acute  if  in 
addition,  broadband  filters  are  used  in  the  analysis.  For  for¬ 
ward  flight,  as  with  hover,  the  narrower  the  filter,  the  better 
the  resolution  to  be  obtained  from  a  given  record. 
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EVALUATION  OF  VORTEX  GENERATED  NOISE 


GENERAL  DISCUSSION 

Since  the  shedding  of  Harman-type  vortices  off  an  airfoil  causes 
variation  of  lift  and  drag  on  the  airfoil  surfaces,  the  result¬ 
ing  oscillatory  force  on  the  airfoil  will  create  noise.  A  vari¬ 
ety  of  empirical  techniques  exist  for  estimating  the  magnitude 
of  noise  caused  by  shed  vortices  from  rotating  wings.  (See 
References  2,  3  and  9.)  The  analysis  carried  out  under  this 
program  determined  the  magnitude  of  these  oscillatory  forces 
from  the  measured  pressure-time  history  of  the  recorded  sound. 

The  RASA  electronically  digitized  pressure-time  histories  were 
separated  into  broadband  and  discrete  components  using  spectral 
analysis  techniques.  The  broadband  pressure-time  history  was 
then  used  in  a  computer  program  which  determined  the  oscillatory 
force  on  the  airfoil  that  created  the  measured  signal.  Because 
of  the  importance  of  the  value  of  the  Strouhal  number  used  in 
the  calculations  and  because  the  Strouhal  number  has  not  been 
accurately  measured,  particularly  for  airfoil  sections,  the 
fitting  technique  used  to  determine  the  oscillatory  forces  was 
also  used  to  determine  the  "best"  Strouhal  number. 

The  data  used  in  this  program  were  only  sufficient  to  determine 
the  contribution  of  each  airfoil  section  to  the  pressure-time 
history.  The  data  were  not  sufficient  to  evaluate  the  relative 
contribution  of  both  lift  and  drag.  Since  it  has  been  estimated 
that  the  oscillatory  drag  forces  are  very  small  compared  to  the 
oscillatory  lift  forces  as  regards  the  acoustic  signal,  the  con¬ 
tribution  from  each  station  was  assumed  to  result  only  from  the 
oscillatory  lift. 


The  Vortex  Shedding  Model 


The  vortex  street  phenomenon  requires  a  separated  wake  caused 
by  the  viscous  effects  in  the  boundary  layer.  While  airfoil 
design  minimized  this  separation  effect,  it  did  not  eliminate 
it.  On  a  physical  basis,  the  wake  formed  from  a  separated 
boundary  layer  will  be  influenced  by  the  wake  width  at  the 
separation  point  and  the  flow  conditions.  The  Strouhal  number 
which  correlates  the  frequency  of  shedding  for  various  bodies 
has  been  defined  as 


(14) 


where  =  Strouhal  number 

u>  *  vortex  street  shedding  frequency  (Hz) 

U  *  flow  velocity  (ft/sec) 
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d  ■  projected  dimension  of  the  body  perpendicular  to 
the  flow  velocity  (ft) 

For  cylinders,  the  Strouhal  number  changes  with  Reynolds  number 
and  is  not  known  for  airfoils  or  cylinders  «.'t  large  Reynolds 
number.  The  Strouhal  number  defined  above  Is  not  a  universal 
constant,  but  is  dependent  on  body  geometry  as  well  as  Reynolds 
number.  For  this  reason,  Roshko  (13)  defined  a  "universal 
Strouhal  mnaber”,  S^,  which  is  not  a  function  of  body  geometry: 


where  SI  -  Roshko 's  "universal  Strouhal  number",  with 
numerical  value  of  about  0.17 

u  -  vortex  shedding  frequency  (Hz) 

U.  -  velocity  just  outside  the  boundary  layer  at  the 
separation  point,  and 

h  -  distance  between  rows  of  vortices  after  the  vortex 
street  is  formed. 


For  thin  streamlined  bodies,  »  U;  data  from  NACA  65  (216)-222 

airfoil  at  angles  of  attack  of  8.1*  and  10.1*  where  separation 
occurred  near  the  maximum  thickness  indicate*;  chat  U.  *  1.14U 

D 

(see  Reference  14).  Therefore,  may  be  considered  bounded  by 
U  and  (1.14)U  for  the  present  analysis. 

The  wake  thickness  h  varies  as  the  projected  dimension  of  the 
body  d  perpendicular  to  the  local  resultant  velocity.  Refer¬ 
ence  3  used  h  *  d  while  Reference  10  used  h  -  (1.54)d.  There¬ 
fore,  Equations  (14)  and  (15)  may  be  rewritten  as 

^11 

s.  -  s;  -  s;  a  (16) 

t  t  Ad  t 

where  An,  Ad,  and  A  are  constants  with  values  defined  by 

1.00  -  AD  -  1.14 
1.00  -  Ad  -  1.54 
0.65  -  A  -  1.14 

Since  S£  -  0.17,  the  Strouhal  number  for  airfoils  can  be  con¬ 
sidered  to  vary  as  follows: 

0.11  -  Sfc  -  0.2  (17) 
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This  range  is  only  approximate,  but  will  be  used  as  a  guide 
until  tests  can  determine  the  existence  of  Strouhal-type  shed¬ 
ding  from  airfoils  and  the  proper  range  and  dependencies  of  the 
strouhal  numbor  on  the  controlling  parameters. 

In  the  present  analysis,  the  Strouhal  number  for  airfoil  vor¬ 
tex  shedding  will  be  allowed  to  vary  within  the  range  0.1  -  Sfc 

-  0.3,  which  excoeds  the  limits  defined  by  Equation  (17).  The 
Strouhal  number  will  be  considered  constant,  however,  along  the 
blade  for  each  calculation;  that  is,  the  Strouhal  number  will 
not  bo  allowed  to  vary  with  blade  velocity  and/or  angle  of 
attack.  It  is  anticipated  that  the  Strouhal  number  is  depen¬ 
dent  on  theso  parameters,  but  in  liou  of  this  knowledge,  the 
approach  used  is  believed  to  be  reasonable.  Thus,  the  vortex 
streot  shodding  frequency  is  assumed  to  be 


(18) 


whoro  u  -  vortex  shedding  frequency  (Hz) 

U  ■  resultant  volocity  (ft/sec) 

d  -  projoctod  dimension  of  the  body  perpendicular  to 
the  resultant  velocity  (ft) 

St  -  Strouhal  number:  0.1  -  S,  -  0.3 


Effects  of  vortox  shodding  on  the  two-dimensional  lift  and  drag 
forces  acting  on  tho  airfoil  can  be  investigated  by  considering 
the  circulation  around  the  appropriate  airfoil  section.  Assum¬ 
ing  rogular  sheddinq  of  vortices  of  alternating  sign,  an  oscil¬ 
latory  circulation  nay  be  considered  to  bo  superimposed  upon 
the  steady  circulation.  The  corresponding  lift  and  drag  forces 
acting  on  the  airfoil  section  will  then  oscillate  about  a  mean 
value  of  lift  and  drag.  The  oscillatory  lift  and  drag  due  to 
vortox  shedding  can  therefore  be  formulated  as 


(Vij(t)  "  <KLv)ij{|1/2pUij]CiAri/  8in  (wijt  4  *ij} 

(19) 

(Dv)ij(t)  "  ^  kdv  >  ij{(1//2puij)ciAri}  sin  Wj*  +  *ij> 


where  (L  ){.(t)  -  oscillatory  lift  at  the  i —  radial  station 
v  1 J  th 

and  tho  j —  azimuthal  station  duo  to  vortex 
shedding  (lb) 

(D  )  >  .  (t)  ■  oscillatory  drag  at  the  i—  radial  station 
v  1 3  th 

and  the  j —  azimuthal  station  due  to  vortex 
shedding  (lb) 
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(K*  *  nondimens ional  force  constant  associated 

13  th 

•  with  oscillatory  lift  at  the  i—  radial 

station  and  the  j  —  azimuthal  station 

■  nondimens  ional  force  constant  associated 

th 

v/ith  oscillatory  drag  at  the  i —  radial 
station  and  the  j  —  azimuthal  station 


p 


“ij 

♦u 


density  of  air  at  the  test  conditions 
lb-sec? 


-  0.002378 


ft'* 


relative  velocity  with  respect  to  the  i — 

radial  station  and  j—  azimuthal  station  of 
the  airfoil  (ft/sec) 

Hi 

chord  length  of  the  i—  radial  station  (ft) 
station  width  of  the  i—  radial  station  (ft) 

frequency  of  vortex  shedding  at  the  i— 

radial  station  and  j—  azimuthal  station 
(rad/sec) 

phase  of  the  vortex  forces  at  the  i— 

til 

radial  station  and  the  j —  azimuthal  station 
(rad) 


The  Cartesian  coordinate  system  is  coincident  with  the  array  of 
microphones  at  which  data  was  recorded.  A  right-handed  coor¬ 
dinate  system  is  used  in  which  the  positive  X-axis  is  in  the 
direction  of  flight,  the  positive  Y-axis  is  to  the  left  of  the 
flight  path,  the  Z-axis  is  positive  upward,  and  the  origin  of 
coordinates  is  at  the  intersection  of  the  microphone  arrays 
(see  Figure  4).  A  local  coordindate  system  (x,y,z)  with  origin 
at  the  helicopter  main  rotor  hub  moves  with  the  helicopter;  see 
Figure  5.  The  positive  x-axis,  however,  remains  in  the  direction 
of  horizontal  flight  with  the  z-axis  p  sitivo  upward. 

Tip-path  coordinates  arc  used  to  locate  a  blade  element  in  the 
x,y,z  coordinate  system.  For  simplicity,  the  hub  is  assumed  to  be 
parallel  to  tho  tip-path  plane.  Rigid  blade  coning  is  accounted 
for  in  locating  a  blade  clement,  hut  the  analysis  neglects  lead- 
lag  motion  and  flexible  blade  bending  in  both  flapwisc  and 
chordwise  directions.  A  blade  element  is  therefore  located 
in  x,y,z  coordinates  as  shown  in  Figure  5  in  which 
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■  -  e  co*  aT  cos  -  (r^  -  e)  co*  60  008  aT  008  ♦j 

■f  (r^  -  e)  sin  Bo  sin  <*T 


yij  “  “  e  *in  *j  “  ^ri  "  008  ®0  8in  *j  f  (20) 

Zij  “  e  sin  aT  008  +  ^ri  “  008  e°  s*n  “t  003 

♦  (ri  -  e)  sin  B0  cos 

where  r ,  *  redial  distance  of  the  i—  blade  element  from  the 
axis  of  rotation  (ft) 

otg,  ■  rotor  angle  of  attack  (positive  for  forward  tixt  of 
tip  plane)  (rad) 

B o  J  coning  angle  relative  to  tip-path  plane  (rad) 

e  ■  distance  from  flapping  hinge  to  axis  of  rotation 
(ft) 

i> .  »  blade  azimuth  angle,  measured  from  negative  x-axis 
3  in  direction  of  rotation  around  tip-path  plane  to 

the  j—  azimuthal  station 
The  location  of  any  source  point  at  time  t  is 

Xij(t)  -  XH(t)  +  x...(t)  (21) 

where  $H(t)  is  the  hub  location  at  time  t. 

The  noise  sources  treated  are  those  generated  from  helicopter 
rotors,  it  is  assumed  that  helicopter  rotor  noise  is  produced 
by  variable  forces  on  the  rotor  blades,  not  fluid  mass  dis¬ 
placement  by  the  blades  nor  flow  turbulence.  It  is  therefore 
reasonable  to  use  dipoles  to  represent  these  blade  forces 
(References  5  and  16).  The  pressure  received  at  an  observer's 
location  from  an  array  of  moving  dipoles  can  be  cast  in  the 
following  form  (see  Reference  5) : 
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V.  . 

ID 


*  velocity  of  blade  element  i,j  (ft/sec) 


f.  .  = 

ID 


M.  . 

ID 


3E 


The  quantities  in  Equation  (?2)  inside  the  brackets  should  be 
evaluated  at  a  retarded  time  prior  to  t  since  the  signal  was 
actually  emitted  at  an  earlier  time  due  to  the  finite  speed  oi 
sound.  That  is,  a  signal  received  at  an  observer's  location 
that  was  emitted  from  a  source  at  time  t'  required  time  At=R/c 
to  reach  an  observer  (where  R  is  the  distance  from  source  to 
observer,  c  is  the  speed  of  sound);  t'  +  At  =  t,  the  real  time 
at  which  the  signal  is  received  by  the  observer. 


If  e  ,  e  ,  and  e  are  unit  vectors  along  the  x,  y,  and  z-axis 
x  y  z 

as  well  as  the  X,  V,  and  Z-axis,  then  the  velocity  of  blade 
element  i,j  is 


v 


ID 


=  (ft 


XX..)  - 


<Vij 


(23) 


where  x..  =  x  ■  •  e  +y..  e  +z..  e 
ID  ID  x  *13  y  id  z 


ft  =  ft I (sin  aT)  ex  +  (cos  aT)  e  ] 


(V  ) . .  =  flow  velocity  due  to  dynamic  response  of  the 
u  13  blade 

ft  =  rotor  rotational  speed 


Similarly, 


v.  .  =  ft  X  (ft  X  *.  .) 

1  J  X  J 


<Vij 


(24) 


since  the  helicopter  is  assumed  to  be  in  steady  rectilinear 
flight  with  constant  rotor  RPM. 


The  inflow  angle  at  blade  element  i , j  is  defined  by 


tan’‘ 

T  1  j 


(25) 
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nri  +  VH  008  °T  s*n  =  con'Ponent  of  resultant 

flow  velocity  perpendicular  to  blade  span  axis 
and  axis  of  rotation 

(VI ) i j  cos  80  +  VH  sin  c*T  cos  B0  +  VH  cos  aT 
times  cos  sin  So  *  component  of  resultant  flow 
velocity  perpendicular  to  blade  span  axis  and 
UT 

induced  inflow  velocity  at  rotor  blade  element 

i/ j 


Following  the  usual  convention,  lift  is  taken  perpendicular  to 
the  resultant  velocity  and  drag  is  taken  parallel  to  the 

resultant  velocity.  Therefore,  the  blade  element  force  com¬ 
ponents  are  defined  by 

?ij  *®x=  {Fij)x=  1  <V  ijcos*i  j-  (V  i jsin*ij 1  (cose0sinaT 

+sinB  ocosaTcos^ j ) 

-[  (Lv}ijsin^ij+  (IVijCOS<*ij]  (COSaTsin'<'j> 

fij-ey=(Fij)y”I(Lv)ijoos*ij-(Dv)ijsln*ij)  (sinBosin^) 

+  t(Lv)ijsii^ij+(Dv)ijcos4»ij]  ( cos i|<  j  ) 

?ij  *ez=  (Fij }  z=  1  (V  ijcos*i  j-  (DV)  i  jsin*i  j]  (cosS0cosaT 

-sinSosinaTcos^ j ) 

+  1  (Lv)i.sin*i:.+  (Dv)ijcos4.i.]  (sinc^sin^) 


where  (UT) ^ 


(V« 


27 


Taking  the  derivatives  with  respect  to  time,  Equation  (26) 
becomes 


(Fij)x-I (Lv)ijC084ij_(^v)ijsin4ij1 <cos*osin3T 

♦  s  in  6  o  co  8 “  fco 8  *  j ) 

-l  (Ly)  ijsin^ij,*’(Dv)  ijcoaiij]  (cosaTsini|(  j  ) 

+  1“  (Ly)  ijSin^ij  -  (Dv)  A  jCOSi)ii  j )  (cosB0sinaT 

+  sin8  ocosa^cosij-  j ) 

(Lv)ijCOS^ij"(Dv)ijsin^ijI  (^ij}  lcos*Tsin* j > 

+  1  (Ly)  ..cos^..-  (Dy)  i^sin^ij  ]  (-nsinB0cosaTsimi<.. ) 

-[  (Lv)  i jSin<^i  j+  (Dy)  iJcosiijJ  (!icosaTcos*. ) 


(^ij)y*  I  (Ly)  ijCOS4ii;.-  (Dy)  ijSin^i  j  ]  (sinB0sin^) 

+  [  (Ly)  ijSin4>ij  +  (Dy)  ijCOSiJ>ij  ]  (cos^^) 

+I-(Lv)ijsin*ij'(Dv)  ijCos^ij]  (^j)  (sinB0sinij'j) 
+  [  (Ly)  ijCOs$i  ^  (Dy)  ijSin<j>ij  ]  (-^j)  (cos*.) 

+  [  (Lv)  ijCQS^j- (Dy)  J  (s]sinB0cos^j ) 

+  [  (Ly)ijSin4>ij+(Dy)ijCOs<^ij]  (-nsimj  j) 

i  jCOs^j  -  (Dy)  ijSin^j  ]  (cosB0cosaT 
-sinBosinaTcosi{' j ) 

+  1  (Ly)  ijSin*i^+  (Dy)  i  jCOs$  i  j  ]  (sinaTsin^.j) 

+  I_(Lv)ijsin<^ij"{Dv)ijcos^ij]  (*ij}  (cosSoCOSotT 

-sinB  0sinaTcos* j ) 

+  1  (Lv)  ijCOS^i;.- (Dy)  ijSin^ij]  (*i;.)  (sinaTsini(j  j ) 

+  1  (Ly)  ^cos^j-  (Dy)  ^sin*^]  (nsinBosinaTsinii» j ) 
+  1  (Ly)  ijSin4>ij+  (Dy)  ijCOs^ij]  (ttsinaTcos* .. ) 

where  (Lv) ^ ^  and  (Dy)ij  are  defined  by  Equation  ( 19)  • 


(27) 
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With  terras  comprising  Equation  (22)  formulated,  the  pressure  at 
an  observer's  location  may  be  specified  once  the  appropriate 
sources  have  been  located  in  retarded  time.  The  acoustic  sig¬ 
nature  of  the  rotor  is  assumed  to  emanate  from  a  series  of 
dipoles  located  at  blade  stations.  The  sound  reaching  a  micro¬ 
phone  at  some  instant  in  tirae,  t,  originates  from  each  of  these 
stations  at  some  earlier  time,  t'.  The  time  that  the  sound 
originated  from  each  station  (retarded  time)  is  a  function  of 
the  position  of  each  station,  and  its  motion  and  is,  in  general, 
different  for  each  station. 

To  determine  the  retarded  time  of  each  station  for  a  time,  t, 
the  following  techniques  were  used.  The  position  <^f  a  station 
on  the  blade  is  known  as  a  function  of  time.  Let  (t)  de¬ 
scribe  the  motion  of  this  station.  Let  X0  be  the  position  of 
the  microphone.  Then  a  sound  wave  originating  from  this  station 
at  some  retarded  time  t'  arrives  at  the  microphone  at  time  t. 

The  time  lapse  of  the  sound  from  inception  to  reception  is  t-t'. 
This  delay  corresponds  to  the  time  required  for  the  sound  wave 
(traveling  at  the  speed  of  sound  c)  to  travel  the  distance 

R  -  |*o  "  *i;)  (f)  |. 


That  is,  jx0  -  ^ ( t ' )  |  =  c (t  -  t') 


(28) 


Equation  (28)  is,  in  general,  transcendental  and  is  solved  using 
an  iterative  technique.  An  initial  estimate  is  made  for  t'. 

This  estimate  t'0  is  used  in  the  left  side  of  Equation  (28)  which 

is  solved  for  the  first  iterative  value,  t'. 

t;  =  t  -  | X0  -  X.j  (t;)  l/c  (29) 

This  new  value  tj  is  then  used  in  the  right  side  of  Equation 

(29),  giving  a  second  iterative  value  t!,  and  so  on  until  the 

solution  converges.  The  retarded  time  calculation,  therefore, 
locates  the  sources  in  space-time  that  contribute  to  the  pres¬ 
sure  at  time  t.  The  retarded  time  calculation  procedure  has 
been  programmed  for  computer  use.  The  program  developed  locates 
for  an  observer's  location  at  time  t,  the  contributing  retarded 
time  position  of  each  radial  blade  station.  The  results  of  a 
typical  retarded  time  calculation  are  shown  in  graphical  form  in 
Figure  32.  The  rotor  disk  has  been  broken  up  into  40  azimuthal 
sections  and  10  radial  sections,  and  the  azimuthal  spacing  in 
retarded  time  is  not  equal  as  it  would  be  in  true  time.  The 
difference  in  the  spacing  between  the  advancing  blade  (observer 
on  the  negative  y-axis)  and  the  retreating  blade  is  very 
obvious  in  this  graphical  presentation. 
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HOVER  ANALYSIS 


For  this  analysis ,  the  UH-lB  helicopter  was  in  the  hover  con¬ 
figuration  100  ft  over  the  microphone  data  cross  with  helicopter 
orientation  along  the  flight  path.  The  data  received  were  re¬ 
corded  at  microphone  4,  which  was  200  ft  to  the  right  of  the 
flight  path. 

Determination  of  Angles  of  Attack,  Resultant  Velocities, 

Shedding  Frequencies 

The  Strouhal  shedding  frequency  of  the  UH-lB  airfoil  (NACA 
0012)  was  assumed  to  be  proportional  to  the  flow  velocity  and 
inversely  proportional  to  the  separation  thickness  of  the  flow. 
This  flow  separation  has  been  assumed  to  be  equal  to  the  pro¬ 
jected  dimension/  d,  of  the  airfoil  perpendicular  to  the  re¬ 
sultant  velocity  U  (see  Equations  (14)  through  (18)).  The 
dependence  of  the  projected  dimension  as  a  function  of  angle  of 
attack  for  the  NACA  0012  blade  is  shown  in  Figure  33.  For  use 
in  computer  programs,  this  curve  was  fit  in  the  three  regions 
indicated  in  Figure  33.  The  functional  relationships  used  in 
each  of  the  three  regions  are: 


where 


Region 


a  Range 


Projected  di 


I 

II 

III 

a  .  is 

l 

C .  is 
a. 


the 

the 


0° 

-  7.5° 

(1  +  0.00623a?)  (0.12 

7.5° 

-  20° 

(0.4  +  0.125ai)  (.12Ci) 

O 

O 

CN 

-  90° 

sin  (ai) 

angle  of  attack  in  radians 
chord  length 


(30) 


For  the  hover  cases  analyzed, the  blade  element  velocities  and 
angles  of  attack  were  calculated  using  the  formulations  of 
Reference  15.  The  wind  velocity  was  small  and  negligible  com¬ 
pared  with  the  blade  velocity  and  therefore  was  assumed  to  be 
zero.  The  inflow  velocity  perpendicular  to  the  blade-element 
direction  of  motion  (Vj)  was 

_  r  i  l'/2l 


<vl>i  - 


n.  c.  nc0 

.  b  l  i-a 


16T 


-  1  + 


1  + 


2»ri (G0+6Ti) 

.  QC  „  1 

l  &a 

6  ti  J 


16i 


(31) 


30 


where  N. 


Ci 


'la 


BTi 


number  of  blades  =  2 

chord  of  i—  station  =  1.75  ft 

main  rotor  rotational  speed  =  33.929  rad/sec 

slope  of  the  curve  of  section  lift  coefficient 
against  angle  of  attack  -  5.73  for  an  NACA  0012 

blade  collective  pitch  angle  =  0.236  rad 

blade  twist  =  -  0.175  r^/rT 

main  rotor  tip  radius  =  22  ft 


The  angle  of  attack  at  the  i—  radial  station  for  hover  is 
defined  by 


=  0n  +  B 


Ti 


-  <J>. 


(32) 


The  flow  velocity  components  transverse  and  parallel  to  the 
blade-element  direction  of  motion  are  respectively  (see  Equation 
(25)) 


where  the  coning 
of  attack  aT  =  0 

The  inflow  angle 


(UT)i  =  nri  and 

(Up)  i  =  (VI)i  cos  B o  (33) 

angle  Bo  =  0*044  radians  and  the  rotor  angle 
for  hover. 

i> .  is  then 


=  tan 


-l 


f<Vi 


“i 


COS  0  o 


(34) 


and  the  magnitude  of  the  total  resultant  flow  velocity  is 

Ui  =  [n2r?  +  (V?  cos 2 B o ]  1//2 


(35) 


Because  the  helicopter  is  in  hover  all  of  these  quantities  are 
assumed  to  have  no  azimuthal  variations. 
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Numerical  Technique 


Measured  pressures  can  now  be  related  to  the  oscillatory  force 
constants.  The  objective  is  to  solve  for  the  unknown  oscil¬ 
latory  force  constants  in  terms  of  the  measured  pressure-time 
history.  For  the  hover  configuration,  the  values  of  the  oscil¬ 
latory  forces  are  assumed  constant  about  the  rotor  disk,  since 
angle  of  attack  and  velocity  are  similarly  assumed  to  be  con¬ 
stant.  The  geometric  parameters,  however,  change  considerably 
with  respect  to  an  observer,  so  it  is  desirable  to  segment  the 
time  for  each  blade  passage  into  smaller  time  intervals.  If 
the  time  intervals  are  small  enough,  the  geometric  parameters 
are  essentially  constant  during  the  interval. 

If  the  geometric  parameters  can  be  assumed  constant  in  a  small 
interval  of  time,  then  the  predictive  equations  can  be  linear¬ 
ized.  When  possible,  such  a  linearization  of  a  set  of  equations 
greatly  reduces  the  numerical  difficulties  of  their  solution. 

The  time  interval  chosen  constrains  the  frequency  range  and 
number  of  oscillatory  stations  per  blade  much  in  the  same  way 
that  length  of  record  is  related  to  bandwidth.  The  largest  time 
interval  that  was  felt  to  be  consistent  with  maintaining  constant 
geometric  parameters  within  that  interval  was  1/200  sec.  For 
example,  if  the  time  for  one  blade  passage  is  divided  into  20 
equal  time  intervals  (the  blade  passage  frequency  is  10.8  Hz), 
each  time  interval  is  approximately  1/200  sec.  Thus,  phenomena 
which  change  with  frequency  less  than  200  Hz  cannot  be  considered. 
This  cutoff  is  a  judicious  compromise  between  frequency  range 
and  blade  element  station  size.  That  is,  if  a  smaller  time  in¬ 
terval  were  chosen,  the  area  swept  by  the  blade  (the  station  size) 
would  be  smaller  but  the  lower  frequency  cutoff  would  increase. 
Conversely,  if  a  larger  time  interval  were  chosen  the  lower  fre¬ 
quency  cutoff  would  decrease,  but  the  area  swept  by  the  blade 
would  increase  and  tend  to  invalidate  this  discrete  element 
analysis . 


The  Strouhal  frequency  formulation  is  such  that  the  shedding 
frequency  increases  with  radius.  Because  of  the  Doppler  effect, 
the  observed  frequencies  are  increased  further  on  the  advancing 
blade  while  decreased  on  the  retreating  blade.  This  is  coupled 
to  the  Mach  number  effect  in  Equation  (22)  which  increases  the 
pressure  amplitude  of  the  advancing  blade  stations  through  the 
terms  1/1-M^  and  decreases  the  pressure  amplitude  of  the  re¬ 
treating  blade  stations  (MR  is  negative  for  a  retreating  blade, 


positive  for  an  advancing  blade) . 
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The  larger  frequency  spread  on  the  advancing  blade  allows  a 
larger  number  of  blade  elements  with  observer  frequency  separa¬ 
tions  above  ^200  Hz.  On  the  retreating  blade  the  Doppler  effect 
works  in  reverse  and  decreases  the  observed  frequencies.  The 
net  result  is  that  for  the  observer  frequencies  above  500  Hz, 
the  section  of  the  airfoil  from  midspan  to  the  tip  of  the  ad¬ 
vancing  UH-lB  blade  contributes,  while  below  500  Hz  the  inboard 
half  of  the  advancing  blade  and  the  whole  retreating  blade 
contribute. 

Solution  Technique 

In  order  to  solve  for  the  magnitude  of  the  oscillatory  forces 
on  the  helicopter  blade  several  assumptions  have  been  made: 


1.  The  oscillatory  forces  on  an  airfoil  occur  at  the 
frequencies  associated  with  Strouhal  shedding. 

2.  The  Strouhal  number  is  assumed  to  be  independent  of 
Reynolds  number. 

3.  The  oscillatory  forces  are  sinusoidal  and  have  com¬ 
ponents  only  in  the  lift  and  drag  direction. 

4.  The  oscillatory  forces  on  the  blade  can  be  represented 
by  oscillatory  dipoles  acting  at  the  center  of  ten 
spanwise  stcitions. 

5.  For  each  small  increment  of  time  ("  0.005  second),  the 
aerodynamic  and  geometric  parameters  at  each  of  the 
ten  spanwise  regions  remain  essentially  constant. 

6.  In  the  time  intervals  when  one  blade  advances  toward 
the  observer  while  the  opposite  blade  retreats,  the 
noise  at  the  observer  for  frequencies  above  500  Hz 
is  assumed  to  originate  from  the  advancing  blade 
above . 


Using  these  assumptions, the  pressure-time  history  at  an  ob¬ 
server  for  a  blade  in  the  advancing  region  is 


P  (t)  = 


i  o 

l 

i=i 


(Kr„).[(AT).  sin  (o).i.+  >|>.)  +  ( B^)  i  cos  (w.i4+i|».)] 


Lv'  i 


L'  l 


L  1 


i  i 


(KDv)iI(AD)i  sin  WW  +  (Vi  COS  (uJiV  V1 


(36) 
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where  (K_  ) .  is  the  magnitude  of  the  oscillatory  lift  coeffi- 
cient  at  the  i—  spanwise  station 

(KDv) i  is  the  magnitude  of  the  oscillatory  drag  coeffi¬ 
cient  at  the  i—  spanwise  station 

is  the  retarded  time  for  station  i 

u).  is  the  Strouhal  shedding  frequency  of  the  i— 
spanwise  station 

is  the  phase  of  the  shedding 

(A_  )  .  ,  (BT  ) .  ,  (A_)  .  and  (B_).  are  the  geometric  coefficients 

L  1  L  1  U  1  U  X 

predicted  by  the  moving,  oscillatory  dipole  theory 
where 
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D  1  «.cK|  !■«,)'  L  i  D1J 


sin  *  sin  4  ♦  cos  *  sin  6  cos  4 
cos  4  sin  4  ♦  sin  4  sin  0  cos  4 
cos  8  cos  4 


?L  *  fl 


sin  4  cos  4  -  cos  4  sin  8  sin  4 

cos  4  cos  4  -  sin  i>  sin  8  sin  4 

cos  8  sin  4 

cos  4  sin  4  -  sin  4  sin  8  cos  4 

sin  4  sin  4  +  cos  *  sin  8  cos  4 

0 


?D  "  n 


-  cos  *  cos  *  +  sin  *  sin  8  sin  4 

-  sin  *  cos  4  -  cos  t  sin  8  sin  4 

0 


(37) 


This  equation  can  be  further  simplified  to 
10 


p(t)  -  l 
i»i 

("V 

^  sin 

(u> 

i’i  * 

4±)  + 

where  (Kg)  ^  ■ 
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<Vi 
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<KC>1  - 

<KL>i 

(BL'i 

■f 

'Vi 

(BD)  i 

Vi 


Vi 


V 

(38) 


The  assumption  that  the  aerodynamic  and  geometric  parameters  re¬ 
main  constant  in  a  0.005-second  time  slice  is  now  applied.  The 
term  up  is  replaced  by  w't,  where  w'  is  the  Doppler  shifted 
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frequency  at  the  observer's  location.  Equation  (38)  then  be¬ 
comes  for  each  blade 

1  o 


P(t)  *  l  <Vi  Bin  («it  ♦  *i)  ♦  (Kc)i  cos  («J 


t  ♦ 


♦i» 


(39) 


The  exact  time  interval  chosen  to  solve  for  the  oscillatory 
forces  was  0.0046  second.  The  rate  at  which  the  data  was 
electronically  digitized  was  11.3  khz.  Therefore,  there  aro 
52  measured  values  of  the  recorded  pressure  in  each  0. 0046- 
second  time  interval.  Linearizing  Equation  (39)  by  removing 
tho  arbitrary  phase  given 


r— 

P(tt)  -  f  jX'Ji  sin  («;tt)  ♦  (K^)i  con  (*[tt)J 


(40) 


t  “  1 , . . .  , 52 


where  /(K^)  *  ♦  (K^)J  .  /ks)J  ♦  (Kc)  |  -  K£ 


and  is  the  total  magnitude  of  the  oscillatory  force  constants. 
Since  there  are  20  unknowns  in  Equation  (40),  ton(Kg)1  and  ten 
(Kc)  £,  and  52  equations,  the  problem  is  overdeterminod .  The  solu¬ 
tion  was  determined  using  the  standard  least  squares  techniques 
for  a  set  of  linear  equations. 

Obtaining  Experimental  Main  Rotor  Vortex  Noise 

The  acoustical  pressure-time  history  recorded  at  the  microphone 
position  is  assumed  to  be  a  linear  superposition  of  many  sources 
of  noiso.  Among  theso  sources  are: 

1.  Engino  noiso 

2.  Engine  exhaust  noiso 

3.  Main  and  tail  rotor  gearbox  noise 

4.  Tail  rotor  rotational  noise 

5.  Main  rotor  rotational  noise 

6.  Main  rotor  vortex  noise 

7.  Tail  rotor  vortex  noise 

Since  the  analyses  that  were  performed  herein  are  associated 
with  only  main  rotor  vortex  noise,  item  (6),  all  the  other 
sources  of  noise  in  the  measured  pressure-time  histories  had 
to  be  removed  in  order  to  obtain  the  desired  pressure-time 
history. 
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In  the  model  developed  in  this  contract,  the  noise  associated 
with  rotor  blade  vortex  shedding  is  assumed  to  be  generated  at 
discrete  frequencies  defined  by  a  Strouhal  formulation.  These 
shedding  frequencies  are  determined  by  the  airfoil's  chord, 
thickness,  angle  of  attack  and  velocity.  Shedding  frequencies, 
while  in  part  determined  by  rotor  rotational  speed,  are  not 
integral  harmonics  of  tho  rotor  rotational  frequency.  In  addi¬ 
tion,  since  the  blade  is  moving  and  continually  changing  direc¬ 
tion,  the  Doppler  shift  on  these  frequencies  also  continually 
changes.  This  differs  from  items  1,  2,  3,  4,  and  5  of  the 
listed  sources  of  noise  which  are  observed  at  discrete  frequen¬ 
cies.  This  is  clearly  illustrated  by  Figure  17,  the  spectrum 
of  measured  helicopter  rotor  noise  for  tho  hover  flight  condi¬ 
tion.  The  rotor  rotational  noise  appears  at  discrete  multiples 
of  the  blade  passage  frequency  (10.8  Hz),  and  the  tail  rotational 
noise  occurs  at  multiples  of  its  blade  passage  frequency  (55.1 
Hz).  The  poaks  at  1922  Hz  and  2133  Hz  are  identified  with  the 
tail-rotor  gear  clash  frequencies  of  the  42-deqrcc  and  90-degree 
(second  harmonic)  gears.  The  frequencies  arc  slightly  shifted 
from  the  onos  given  in  Table  II  because  of  the  inability  of  the 
helicopter  to  hold  an  exactly  constant  RPM.  No  other  engine  or 
transmission  noise  is  observed  in  this  data  (see  Table  II). 

This  is  believed  to  be  duo  to  the  generally  higher  frequency 
content  of  those  sources  of  noise  and  the  large  distance  of  the 
microphone  from  tho  helicopter. 

Figure  17  shows  that  the  rotational  and  discrete  noise  peaks 
are  sufficiently  narrow  that  they  can  be  removed  without  signi¬ 
ficantly  changing  the  content  of  tho  underlying  broadband  noise. 
Once  the  discrete  noise  peaks  arc  identified,  their  correspond¬ 
ing  Fourior  coefficients,  as  determined  from  tho  fast  Fourier 
transform  of  the  digitized  data,  are  set  equal  to  zero.  Tho 
broadband  pressure-time  history  is  then  generated  by  using  the 
inverse  Fourier  transform. 

In  order  to  correlate  the  measured  pressure-time  history  with 
the  theory,  the  azimuthal  blade  position  must  be  known.  Tho 
Wallops  Island  test  did  not  rocord  this  important  information. 

Tho  azimuthal  location  of  the  rotor  was  determined  by  matching 
tho  measured  rotational  noi3e  signature  with  a  theoretically 
predicted  rotational  noise  prossuro-time  history. 

Results  o f  the  Theoretical  frnalysis 

Four  hover  cases  wore  analyzed.  Three  of  these  cases  were  from 
position  4  (200  foot  to  tho  right  of  the  helicopter)  and  one 
from  position  6  (700  foot  to  tho  right  of  the  helicopter).  The 
radial  stations  selected  were: 

ri  -  7,  9,  11,  13,  15,  17,  1 R . 5 ,  19.5,  20.5  and  21.5  ft 
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TABLE  II.  GEAR  MESH  FREQUENCIES  OF  UH-1B  DRIVE  TRAIN 
FOR  ENGINE  OUTPUT  SHAFT  SPEED  OF  6600  RPM 


Transmission  Group 

RPM 

Number 

of 

Teeth 

Gear  Mesh 
Frequency 
(Hz) 

Main  Transmission: 

Input  Drive  Shaft  Pinion 

6600 

29 

3190 

1st  Stage  Planetary  Pinions 

3843 

31 

1986 

2nd  Stage  Planetary  Pinions 

1245 

31 

643 

Center  Shaft  Bevel  (Lower) 

3090 

55 

2833 

Tail  Rotor  Output  Shaft  Gear 

4140 

27 

1863 

Generator  Quill 

— 

— 

4510 

Tail  Rotor  Transmission: 

42-Degree  Gearbox  Bevel 

4300 

27 

1935 

90-Degree  Gearbox  Bevel 

4300 

15 

1075 

Engine  Group 

Gearbox : 

1st  Stage  Pinions 

— 

- 

6150 

2nd  Stage  Pinions 

— 

— 

3590 

Accessory  Drive  Gearbox: 

Accessory  Drive  Pinion 

— 

- 

10595 

Accessory  Drive  Gear 

— 

- 

4162 

Oil/Air  Separator  Shaft  Gear 

- 

- 

2973 

Oil  Pump  Drive  Shaft  Gear 

- 

- 

2973 

Tachometer  Drive  Shaft  Gear 

- 

- 

2727 

Fuel  Control  Drive  Shaft  Gear 

• 

2973 

NOTE:  All  blade  passage  frequencies  of 

the  power 

turbine, 

gas  producer,  axial  compressors 

and  centrifugal 

compressor  are  above  1&000 

Hertz 

• 

for  i=l  to  10  respectively.  The  electronically  digitised  re^ 
cords  are  ^1.5  seconds  long.  Therefore,  15  blade  passages  were 
analyzed  for  each  of  the  cases. 

The  set  of  equations  solved  had  more  equations  (52)  than  un¬ 
knowns  (20).  By  using  the  least-squares  method,  the  result¬ 
ing  deviation  is  a  measure  of  the  quality  of  the  solution. 

Since  the  Strouhal  number  is  not  accurately  known  for  airfoils 
at  these  Reynolds  numbers,  a  search  was  performed  to  determine 
the  Strouhal  number  which  gave  the  most  satisfactory  solutions. 
The  optional  Strouhal  number  was  determined  to  be  0.235+.02, 
which  is  within  the  range  for  which  results  have  been  obtained 
by  previous  investigations. 

One  of  the  recordings  obtained  200  ft  to  the  right  was  used 
to  evaluate  the  method  of  removing  discrete  rotational  noise 
sources.  The  portion  of  the  recording  used  was  35  seconds 
into  the  recorded  signal.  The  total  recorded  pressure-time 
history  for  this  portion  of  the  recording  is  shown  in  Figure 
34,  and  the  spectrum  analysis  of  the  signal  is  presented  in 
Figure  35.  In  the  first  of  the  four  cases  (to  be  referred 
to  as  Case  A),  for  all  noise  below  140  Hz,  the  tail  rotor  peaks 
at  165  Hz  and  220  Hz,  and  all  noise  above  1860  Hz  was  re¬ 
moved.  The  spectra  of  the  noise  removed  are  shown  in  Figure 
36,  and  the  remaining  spectrum  that  is  assumed  to  be  associated 
with  vortex  noise  is  shown  in  Figure  37 .  The  pressure-time 
history  corresponding  to  the  discrete  noise  (the  main  and 
tail  rotor  rotational  noise  and  the  tail  rotor  transmission 
noise)  spikes  at  1930  Hz  and  2140  Hz  is  shown  in  Figure  38. 

The  pressure-time  history  of  the  remaining  "vortex"  noise 
for  Case  A  is  shown  in  Figure  39.  Comparison  of  the  "vortex 
noise"  pressure-time  history  (Figure  39)  and  the  discrete 
noise  pressure-time  history  (Figure  38)  with  the  pressure- 
time  history  of  the  total  measured  signal  indicates  that  the 
majority  of  the  high  frequency  content  is  associated  with 
"vortex  noise"  and  not  with  discrete  noise  such  as  main  and 
tail  rotor  rotational  noise.  While  it  is  evident  in  the 
pressure-time  history  of  the  total  noise  signal,  the  increase 
in  the  level  and  amount  of  high  frequency  noise  at  each 
blade  passage  is  very  obvious  and  undi sputable . 

In  Case  B,  all  noise  above  28  20  Hz  war,  removed  (instead  of 
above  1860  Hz  as  in  Case  A),  together  with  the  gear  clash 
spikes  at  around  1910  Hz  and  2140  Hz.  The  spectra  for  the 
removed  noise  and  the  remaining  noise,  as  well  as  the 
pressure-time  history  of  the  romaininq  noise,  are  shown  in 
Figures  40,41  and  42  respectively  Comparing  the  pressure-time 
histories  presented  in  Figures  39  and  42  indicates  that  the 
higher  cutoff  frequency  did  not  notably  change  the  character 
of  the  pressure-time  history  of  the  vortex  noise. 
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In  order  to  determine  how  much  the  results  were  dependent  on 
when  the  data  were  analyzed ,  another  portion  of  the  data  taken 
200  feet  to  the  right  and  26  seconds  into  the  recording  was 
analyzed.  This  condition  is  referred  to  as  Case  C.  The 
electronically  digitized  pressure-time  history  and  spectrum 
for  this  case  are  shown  in  Figures  43  and  17  respectively.  The 
rotational  and  discrete  noise  was  removed  by  setting  the  follow¬ 
ing  coefficients  of  the  Fourier  representation  of  the  digitized 
record  equal  to  zero:  <140  Hz,  165  Hz,  220  Hz,  1930  Hz,  2140  Hz, 
and  >2820  Hz.  Figure  44  presents  the  spectrum  of  the  noise  that 
was  removed,  and  Figure  45  presents  the  spectrum  of  the  vortex 
noise  that  was  analyzed  for  this  case.  The  pressure-time  history 
of  the  rotational  noise  is  shown  in  Figure  46  and  of  the  vortex 
noise  in  Figure  47. 

The  fourth  case  that  was  analyzed  (Case  D)  was  data  measured 
700  feet  to  the  right  of  the  helicopter.  A  portion  of  the 
recording  26  seconds  into  the  record  was  analyzed.  The 
pressure-time  history  and  spectrum  of  this  data  are  shown  in 
Figures  28  and  18  respectively.  In  this  position,  the  gear 
clash  frequencies  occur  at  1930  Hz  and  1990  Hz.  The  peak  at 
1990  Hz  corresponds  to  the  gear  clash  frequency  of  the  first- 
stage  planetary  pinion  (see  Table  II).  The  frequency  com¬ 
ponents  removed  were  <140  Hz,  165  Hz,  220  Hz,  276  Hz,  330  Hz, 

385  Hz,  1930  Hz,  1990  Hz,  and  >2820  Hz.  The  spectra  and 
pressure-time  histories  for  the  noise  removed  and  remaining 
are  presented  in  Figures  49,  50,  51,  and  52.  The  pressure-time 
history  of  the  rotational  noise  is  considerably  different  than 
that  at  the  300  ft  position;  and  the  pressure-time  history  of 
the  vortex  noise,  while  having  the  same  general  characteristics 
of  those  at  200  ft,  has  some  different  features.  The  analyses 
that  were  carried  out  used  15  blade  passages  of  digitized  data. 
Because  of  the  inability  to  separate  the  lift  and  drag  compo¬ 
nents  and  because  the  drag  contribution  is  insignificant  com¬ 
pared  to  the  lift  distribution,  the  (KD)^'s  were  assumed  zero 

in  all  of  the  analyses  that  wore  conducted.  Therefore,  only 
(K.).'s  were  calculated  for  the  15  blade  passages  in  each  case. 

For  each  radial  station,  the  distribution  of  these  15  (K^)^ 

values  multiplied  by  that  station's  dynamic  pressure  is  shown 
in  Figure  53.  The  (K^) ^ (l/2pUp  for  the  ten  radial  stations 

for  each  of  the  four  cases  are  tabulated  in  Table  III  and  plotted 
in  Figure  54.  Although  there  is  some  scatter  in  each  of  the 
distributions,  the  mean  results  are  remarkably  similar  for  all 
four  cases.  The  comparison  of  Case  A  and  Case  B  shows  that  the 
results  are  fairly  insensitive  to  variations  in  the  high  fre¬ 
quency  cutoff  in  the  data,  although  the  magnitudes  of  station 
10  for  Case  A  are  somewhat  suppressed  from  Case  B.  Going  to  a 


40 


41 


different  position  of  the  record  also  does  not  significantly 
alter  the  distributions  (see  distributions  for  Case  C) .  Some 
of  the  distributions  for  position  6  (Case  D) ,  however,  do  differ 
from  the  first  three  cases.  This  is  due,  in  part,  to  the  differ 
ences  in  ground  reflection.  The  magnitudes  of  the  oscillatory 
pressures  are  ^1  to  2  lb/ft2  for  the  tip  1/3  of  the  blade 
(14-22  ft) .  On  tne  inboard  stations  (6-14  ft)  ,  the  magnitude  of 
the  oscillatory  pressure  ranges  from  4  to  30  lb/ft2.  These 
higher  values  associated  with  the  inboard  stations  are  due  in 
part  to  the  fact  that  the  retreating  blade's  range  of  Doppler- 
shifted  shedding  frequencies  is  the  same  as  the  frequency  range 
associated  with  the  inboard  stations  of  the  advancing  blade. 
Therefore,  the  calculated  amplitudes  for  inboard  stations  of 
the  advancing  blade  contain  the  noise  energy  of  the  retreating 
blade.  For  each  of  the  four  cases  that  were  analyzed,  the  cal¬ 
culated  oscillatory  forces  were  used  in  a  theoretical  program 
to  determine  if  the  essential  noise  characteristics  had  been 
retained.  For  each  of  the  cases,  therefore,  a  pressure-time 
history  and  a  corresponding  spectral  analysis  were  predicted 
for  comparison  with  the  experimental  data.  The  predicted 
pressure-time  history  and  the  corresponding  spectrum  for  Case  A 
are  shown  in  Figures  55  and  56  respectively.  When  these  pre¬ 
dicted  results  are  compared  with  the  measured  results  in 
Figures  39  and  37,  it  is  seen  that  the  calculated  oscillatory 
forces  have  retained  all  of  the  essential  characteristics  of 
the  measured  pressure-time  history.  Similar  comparisons  can  be 
made  for  Cases  B,  C  and  D.  For  example,  Figures  57  and  58 
should  be  compared  with  Figures  42  and  41  for  Case  B;  Figures 
59  and  60  should  be  compared  with  Figures  47  and  45  for  Case  C; 
and  Figures  61  and  62  should  be  compared  with  Figures  52  and 
50  for  Case  D. 

Forward  Flight  Analysis 

Calculations  similar  to  those  carried  out  fcr  hover  were  con¬ 
ducted  for  a  high-speed  forward  flight  case  (u=0.28).  The 
helicopter  was  at  100-ft  altitude  and  flew  over  the  data  cross 
in  the  X-direction  (see  Figure  4)  at  208  ft/*ec.  This  problem 
is  considerably  more  difficult  than  hover  for  the  following 
reasons : 

1.  The  helicopter's  position  at  any  instant  of  time  was 
measured  only  to  within  45  ft  (Reference  12). 

2.  The  Doppler  effect  makes  it  more  difficult  to  separate 
the  discrete  rotational  and  gearbox  noise. 

3.  The  rotor  wake  and  blade  loads  are  not  as  easily  or 
accurately  calculated. 
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4.  There  is  an  azimuthal  variation  in  the  aerodynamic 
parameters  such  as  angle  of  attack,  inflow  velocity 
and  flow  velocity. 

To  improve  on  the  uncertainty  of  the  helicopter's  position,  the 
period  of  rotational  noise  as  a  function  of  time  was  matched 
to  a  predicted  Doppler  shift  dependence.  For  these  calculations 
the  helicopter  was  assumed  to  be  in  straight  and  level  flight 
with  altitude  and  velocity  constant  as  given  by  the  radar  mea¬ 
surements.  The  rotor  RPM  was  assumed  to  be  324.  Figure  63 
shows  the  period  of  each  blade  passage  as  a  function  of  time. 

As  expected,  the  period  is  small  where  the  helicopter  is  far 
away  and  approaching,  and  gradually  increases  until  it  is  large 
and  relatively  more  constant  when  the  helicopter  is  far  away 
and  receding.  These  limiting  frequencies  (at  ±2  sec)  from  the 
measured  data  correspond  exactly  to  the  Doppler  shifts  expected 
for  this  flight  configuration.  The  curve  superimposed  on  the 
data  is  the  predicted  period  time  history.  From  this  plot  the 
time  in  the  record  corresponding  to  the  helicopter  crossing  at 
ground  zero  can  be  determined  with  an  error  of  only  one  blade 
passage  or  20  ft,  which  is  considerably  better  than  the  measured 
data. 

The  portion  of  the  record  electronically  digitized  for  the  for¬ 
ward  flight  case  was  chosen  so  that  the  helicopter  was  near  the 
data  cross.  For  this  helicopter  position,  as  in  the  hover  case, 
the  advancing  blade  dominates  the  retreating  blade  noise,  and 
the  Doppler  shift  on  the  advancing  blade  oscillating  dipcles 
enables  the  accurate  separation  of  the  noise  from  the  advancing 
and  retreating  blades.  The  advancing  blade  location  means  the 
blade  is  advancing  toward  the  observer,  and  for  microphone 
position  4,  this  corresponds  approximately  to  c^O0.  The  aero¬ 
dynamic  environment  of  the  blade  for  this  case  is  similar  to 
that  of  hover,  since  the  flow  caused  by  the  helicopter  forward 
velocity  is  only  along  the  blade.  When  the  blade  is  at 
some  frequency  separation  is  expected,  since  the  flow,  and  hence 
the  shedding  frequency,  is  higher  on  that  blade  than  on  the  blade 
at  ij»=270°. 

Figure  64  shows  the  retarded  time  diagram  for  the  two-bladed 
UH-1B  in  forward  flight.  The  values  of  the  wake  inflow  velocity 
were  calculated  using  a  deformed  wake  analysis  d^v  loped  by  RASA 
and  presented  in  Reference  17.  The  /otatioi  *  X  9  was  re" 
moved  using  a  technique  similar  t(  uh.  *  ver.  However, 

because  of  the  broadening  of  the  pc  bdi  spikes  caused  by 

Doppler  effects,  the  periodic  no 1  se  as  ,nved  by  zeroing  out 
the  Fourier  coefficients  below  280  Hz  a  ji/ove  Hz. 

For  the  time  intervals  when  the  b)i(  'e  wus  advancing  toward  the 
observer's  location,  the  reau!’*  c  the.  analyses  were  in  com¬ 
plete  agreement  with  the  hover  data;  therefore  the  hover  results 
also  apply  for  this  forward  flight  case. 
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GENERAL  COMMENTS  AND  SPECIFIC  CONCLUSIONS 


The  technique  of  electronically  digitizing  and  analyzing  mea¬ 
sured  sound  pressure-time  histories  to  create  high  resolution 
spectra,  such  as  shown  in  Figure  17,  has  permitted  the  detailed 
analysis  of  various  rotor  sources.  For  example,  the  fine  fre¬ 
quency  resolution  that  can  be  obtained  by  these  techniques 
allows  separation  of  the  individual  peaks  associated  with  main 
rotor  and  tail  rotor  rotational  noise  and  discrete  sources  such 
as  oil  pumps,  gearboxes,  etc.,  from  the  broadband  noise.  Since 
these  types  of  noise  sources  can  now  be  adequately  separated, 
the  characteristics  of  the  various  noise  sources  can  be  studied 
independently.  The  study  stressed  the  analysis  of  the  broadband 
noise  rather  than  the  rotational  noise,  although  the  rotational 
and  discrete  noise  sources  were  separated  from  the  total  noise 
signature  and  studied  as  regards  the  characteristics  of  their 
pressure-time  history.  Figures  38,  46,  and  51,  for  example, 
show  the  pressure-time  histories  of  the  rotational  and  discrete 
noise  sources  for  the  various  cases  that  were  studied.  While 
the  general  characteristics  of  these  noise  sources  are  the 
same,  there  are  noticeable  differences.  For  example,  comparing 
the  signatures  shown  in  Figures  38  and  46  which  were  recorded 
at  different  times  but  at  the  same  location  during  a  hover 
flight  condition,  it  can  be  seen  that  in  one  case  the  tail  rotor 
signature  was  more  pronounced  than  in  the  other  case.  It  is 
believed  that  this  is  due  to  the  fact  that  the  helicopter  was 
not  holding  a  constant  heading  at  all  times.  Comparing  Figures 
46  and  51,  it  can  be  seen  that  for  a  given  time,  the  effect  of 
microphone  location  can  change  the  character  of  the  pressure¬ 
time  history:  at  700  feet  from  the  aircraft  the  tail  rotor  sig¬ 
nature  was  very  pronounced,  whereas  at  200  feet  it  is  hardly 
distinguishable . 

It  is  apparent  from  this  brief  analysis  that  considerable  care 
must  be  utilized  in  evaluating  the  relative  importance  of  vari¬ 
ous  sources  of  helicopter  noise  based  on  a  given  recording, 
since  many  geometric  and  environmental  parameters  must  be  prop¬ 
erly  considered  in  the  evaluation.  For  any  future  studies  of 
the  type  conducted  herein,  controlled  tests  in  the  hover  condi¬ 
tion  should  be  undertaken  so  that  the  effects  of  blade  loading, 
Mach  number,  etc.,  can  be  adequately  recorded  and  documented. 

The  vortex  pressure-time  histories  (Figures  39,  42,  47,  and  52) 
that  were  generated  from  the  total  experimental  pressure-time 
histories  by  removing  the  rotational  and  other  discrete  noise 
have  the  same  general  characteristics.  The  signals  are  essen¬ 
tially  random,  with  a  modulation  in  amplitude  occurring  every 
blade  passage.  The  high  amplitude  region  has  higher  frequency 
than  the  low  amplitude  regions.  This  feature  is  consistent 
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with  the  Doppler  effect  discussed  previously,  i.e.,  the  advanc¬ 
ing  blade  towards  the  observer  raises  the  frequency  and  magni¬ 
tude  of  any  oscillatory  pressure  amplitudes  recorded  at  an 
observer's  location.  Each  blade  passage,  however,  is  distinctly 
different  in  its  structure.  In  order  to  evaluate  the  effect  of 
this  difference  on  the  sound  signature  of  vortex  noise,  each 
blade  passage  of  the  digitized  vortex  signal  was  converted  to 
an  analog  signal  and  then  repeated  so  that  a  5-second  analog 
record  of  each  particular  blade  passage  could  be  constructed 
on  tape.  When  a  series  of  each  of  these  5-second  records  was 
played  on  a  tape  recorder,  each  blade  passage  sounded  distinctly 
different.  Qualitatively,  this  difference  from  blade  passage 
to  blade  passage  may  be  described  as  a  modulated  signal  of 
varying  frequency  or  "tone".  That  is,  when  listening  to  a  series 
of  blade  passage  recordings,  a  different  frequency  content  may 
be  discerned  in  each  of  the  blade  passages. 

The  vortex  shedding  model  allowed  a  finite  frequency  range  to 
be  "fit"  to  the  experimental  pressure-time  history.  The  rela¬ 
tive  magnitudes  of  each  of  the  radial  station's  oscillatory 
forces  reflected  the  frequency  content  of  the  signal  received  at 
the  observer  for  that  blade  passage.  Thus,  since  the  signal 
varies  from  blade  passage  to  blade  passage,  both  in  frequency 
and  magnitude,  the  oscillatory  pressures  calculated  at  the 
helicopter  blade  had  a  corresponding  variation.  Since  each 
blade  passage  is  different,  each  calculation  performed  results 
in  a  slightly  different  radial  array  of  oscillatory  forces. 

This  effect  is  shown  in  the  histograms  (see  Figure  53) ,  which 
illustrate  how  these  oscillatory  pressures  varied  at  each  radial 
station  over  the  15  blade  passages  for  each  of  the  four  cases 
that  were  analyzed.  Histograms  plot  the  frequency  of  occur¬ 
rences  of  a  given  event.  In  this  case,  the  histograms  record 
the  frequency  of  occurrence  of  the  calculated  oscillatory  lift 
forces  in  a  given  range.  Figure  53  shows  that  the  scatter  of 
results  increases  as  the  radial  station  decreases,  which  relates 
to  the  gradual  degeneracy  of  the  solution  technique  as  the  fre¬ 
quency  of  shedding  on  the  blade  decreases.  The  amount  of 
scatter  of  a  station  corresponds  to  the  variability  of  the 
amplitude  appropriate  to  that  station  contained  in  the  experi¬ 
mental  record.  Thus,  for  Case  C,  the  amplitude  appropriate 
to  station  2  was  changing  significantly  from  blade  passage  to 
blade  passage  while  the  amplitude  appropriate  to  station  6 
changed  little  from  blade  passage  to  blade  passage. 

The  mean  values  for  each  of  these  distributions  were  determined. 
The  mean  values  are  listed  in  Table  III  and  plotted  in  Figure 
54.  The  shedding  frequencies  appropriate  to  the  helicopter 
rotor  reference  system  are  listed,  as  well  as  the  Mach  number 
and  angle  of  attack. 
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The  advancing  blade  has  its  shedding  frequencies  raised  and 
spread  apart,  while  the  retreating  blade  has  its  shedding  fre¬ 
quencies  lowered  and  pushed  together.  The  Doppler  shifted 
frequencies  on  the  advancing  blade  range  from  about  275  Hz  at 
7  ft  radius  to  2100  Hz  at  21.5  ft  radius.  The  exact  amount  of 
the  shift  depends  on  the  observer's  location  relative  to  the 
velocity  of  that  blade  station.  The  Doppler-shifted  frequen¬ 
cies  on  the  retreating  blade  range  from  180  Hz  at  7  ft  radius 
to  about  500  Hz  at  2.15  ft  radius.  Hence,  the  outboard  section 
of  the  advancing  blade  is  associated  with  high-frequency  noise; 
and  the  inboard  section  of  the  advancing  blade,  together  with 
the  retreating  blade,  is  associated  with  low-frequency  noise. 

For  the  observer,  the  radial  stations  at  7  ft  and  9  ft  on  the 
advancing  blade  are  in  the  same  frequency  range  as  the  retreat¬ 
ing  blade  because  of  the  Doppler  effect;  hence,  these  stations 
also  reflect  the  noise  energy  of  the  retreating  blade.  This 
contribution  partially  leads  to  higher  values  of  oscillatory 
pressure  calculated  at  these  radial  stations.  The  noise  energy 
of  the  entire  retreating  blade  has  been  lumped  into  these  two 
stations . 

As  shown  in  Figure  54,  the  decreasing  oscillatory  pressure 
magnitudes  at  larger  blade  radius  reflect  the  gradual  fall-off 
with  frequency  of  the  noise  energy  shown  in  the  spectra  (see 
Figure  17  for  example) .  The  spectra  generated  in  this  analysis 
compare  well  with  the  findings  of  Cox  and  Lynn,  Reference  18, 
in  that  the  major  source  of  audible  vortex  noise  is  concentrated 
in  the  frequency  range  of  200  Hz  to  500  Hz.  The  maxima  and 
minima  in  these  spectra  caused  by  ground  reflection  also  affect 
the  calculated  oscillatory  pressures,  as  the  model  does  not  in¬ 
clude  this  reflection  effect. 

The  mean  oscillatory  pressures  have  been  used  to  create  a 
pressure-time  history  at  an  observer's  station.  These  created 
pressure-time  histo;  ies  and  these  spectra  compare  very  well 
with  the  experimental  pressure-time  histories  and  spectra, 
see  Figures  55,  56,  39  and  37  for  example.  Audio  tapes  created 
from  the  calculated  pressure-time  history  sound  like  the  experi¬ 
mental  tapes  except  for  the  amplitude  variation  from  blade 
passage  to  blade  passage  eliminated  by  using  the  mean  values. 

Since  the  essential  characteristics  of  the  experimental  acoustic 
signal  can  be  duplicated  with  these  mean  oscillatory  forces 
in  a  rotating  blade  frame  of  reference,  it  is  possible  with 
proper  parameterization  of  the  force  constants  to  effectively 
simulate  the  noise  produced  by  helicopter  rotors  for  a  variety 
of  operating  conditions.  Because  of  the  advanced  data  analysis 
procedures  that  have  been  developed,  the  rotational  and  broad¬ 
band  vortex  noise  can  be  separated  and  studied  independently. 

It  is  therefore  possible  to  parameterize  the  vortex  noise  and 
the  rotational  noise  as  functions  of  the  rotor  geometric  and 
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operating  conditions  such  as  number  of  blades,  rotational  speed, 
thrust,  chord,  and  twist.  The  parameterization  of  the  vortex 
noise  and  rotational  noise  can  be  carried  out  by  analysing 
various  rotor  data  from  whirl  tower  tests,  as  was  done  in  this 
program  for  a  hovering  helicopter.  Once  the  parameterization 
has  been  accomplished  for  a  number  of  different  rotor  systems, 
it  should  be  possible  to  predict  the  acoustic  signature  of  any 
rotor  system  given  only  the  geometric  and  operating  conditions 
of  the  helicopter. 

CONCLUSIONS 


The  analyses  carried  out  in  this  investigation  showed  that  ran¬ 
dom  or  "vortex"  noise  is  a  major  source  of  noise  from  a  heli¬ 
copter,  particularly  in  the  hover  condition  For  the  cases 
analyzed,  main  rotor  rotational  noise  was  not  significant  above 
100  Hz,  and  the  major  audible  sources  of  vortex  noise  were  con¬ 
centrated  in  the  frequency  range  of  200  to  500  Hz.  In  their 
respective  frequency  ranges,  the  main  rotor  rotational  noise 
was  about  30  db  above  the  threshold  of  hearing  while  the  vortex 
noise  was  about  45  db  above  a  normal  background  noise.  The 
effect  on  an  acoustic  spectrum  of  the  filter  bandwidth  and 
the  time  length  of  data  was  that  as  the  time  period  decreased 
or  the  filter  bandwidth  increased,  the  resolution  of  the  spec¬ 
trum  decreased. 

The  following  specific  conclusions  were  drawn: 

1.  Ground  reflection  effects  can  significantly  distort 
the  radiated  acoustic  pressure-time  history. 

2.  The  broadband  noise  created  by  a  helicopter  rotor  can 
be  represented  by  Karman-street  vortex  induced  noise. 

3.  The  "vortex"  noise  theory  fit  the  data  best  for  a 
Strouhal  number  of  0.235. 

4.  The  basis  of  a  technique  by  which  the  acoustic  signa¬ 
ture  of  any  rotor  system  can  be  predicted  in  hover  or 
forward  flight  has  been  developed  and  demonstrated. 

5.  A  new  technique  has  been  developed  whereby  jni-fi- 
cantly  improved  high  resolution  acoustic  si.  v  can 
be  generated. 

6.  While  acoustic  spectra  of  helicopter  flybys  cai  u* 
generated,  the  information  that  can  e  obtained  'rom 
such  spectra  is  marginal  because  of  the  effects  of 
Doppler  shift  and  changing  propagation  distance. 
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Figure  1.  Noise  Spectrum  by  Ubiquitous  Analyzer  for  UH-1B  Helicopter  in 
100-ft  Hover;  Recorded  at  200  ft  to  the  Right. 
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Figure  2.  Noise  Spectrum  by  B  &  K  Constant-Percent  Analyzer  for 
UH-1B  Helicopter  in  100-ft  Hover;  Recorded  at  200  ft 
to  the  Right. 


Record  Length 


FREQUENCY  (HZ) 

Figure  3.  Noise  Spectrum  by  Ubiquitous  Analyzer  for  UH-1B  Helicopter 
in  ICO-ft  Hover;  Recorded  at  700  ft  to  the  Right. 
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Figure  12.  Oscillograph  Recording  of  Pressure-Time  History  for 

UH-lB  Approaching  at  115  Knots  Level  Flight;  Recorded 
at  200  ft  to  the  Right. 
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UH-1B  Approaching  at  115  Knots  Level  Flight;  Recorded 
at  200  ft  to  the  Right. 
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Figure  14.  Oscillograph  Recording  of  Pressure-Time  History  for 
UH-1B  in  100  ft-Hover,  Recorded  at  700  ft  to  the 
Right . 
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Figure  16.  Change  in  Character  of  Tail  Rotor  Noise 
During  115-Knot  Level  Flight  Approach. 
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Figure  17.  RASA  Electronically  Digitized  Noise  Spectrum  for 
UH-1B  in  Hover;  Recorded  at  200  ft  to  the  Right; 
Data  Taken  26  Seconds  Into  Record. 
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Figure  18.  RASA  Electronically  Digitized  Spectrum  for  UH-lB 
in  100  ft  Hover;  Recorded  at  700  ft  to  the  Right; 
Data  Taken  26  Seconds  Into  Record. 
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Figure  26.  Analytically  Determined  Noise  Spectrum  at  Position 
for  UH-1B  in  llover;  Length  of  Record  Is  Ten  Blade 
Passages. 
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Figure  31.  RASA  Electronically  Digitized  Noise  Spectrum  at 
Position  4  for  UH-lB  in  Hover. 
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Figure  32.  Retarded  Time  IXtAgram  for  Hover  Showing 
Locations  of  Blade  Radial  and  Azimuthal 
Stations . 
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Figure  38.  RASA  Electronically  Digitized  Pressure-Time  History 

of  Rotational  and  Gear-Clash  Noise  for  UH-1B  in  Hover. 
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Figure  40.  RASA  Electronically  Digitized  Noise  Spectrum  at  Position 
for  UH-1B  in  Hover  (Vortex  Noise  Removed) . 
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Figure  42.  RASA  Electronically  Digitized  Pressure-Time  History  of 
Vortex  Noise  at  Position  4  for  UH-1B  in  Hover. 
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Figure  47.  RASA  Electronically  Digitized  Pressure-Time  History  of 
Vortex  Noise  at  Position  4  for  UH-1B  in  Hover. 
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Figure  48.  RASA  Electronically  Digitized  Pressure-Time  History  of 
Total  Recorded  Signal  at  Position  6  for  UH-1B  in  Hover. 
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Figure  52.  RASA  Electronically  Digitized  Pressure-Time  History  of 
"Vortex"  Noise  at  Position  6  for  UH-1B  in  Hover. 
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Figure  53.  Distribution  of  the  Magnitude  of 
Oscillatory  Lift  Over  Ten  Radial 
Stations  for  Four  Cases. 
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Figure  54.  Magnitudes  of  the  Oscillatory  Lift  Forces  as  a 
Function  of  Radial  Position  for  Each  of  the 
Four  Cases. 
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Figure  55.  Predicted  Vortex  Pressure-Time  History  at  Position 
for  UH-1B  in  Hover. 
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Figure  57.  Predicted  Vortex  Pressure-Time  History  at  Position 
for  UH-iB  in  Hover. 
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UH-lB  in  Hover 
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Figure  61.  Predicted  Vortex  Pressure-Tine  History  at  Position 
for  UH-lB  in  Hover. 
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Figure  63.  Period  of  Rotational  Noise  as  a  Function  of  Time 
at  115  Knots. 
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